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BY ; 
J. KAMPE pe FERIET 
University of Lille, France 
Introduction. In this paper, we consider the two dimensional flow of an incom- 
pressible viscous fluid . 
a) in a bounded domain D, the velocity vanishing at the boundary B of the 
domain; B is assumed to be a regular curve;! 
b) in the entire plane, the velocity vanishing at infinity, this case being regarded 
as a limiting case of a). 
Since the fluid is assumed to be incompressible, there exists a stream function p(x, y, t) 
such that the two components u and v of the velocity are given by 
oy 


u(x, y,t) = a (1) 


oy 
v(x, y, t) = ax F 
x 


According to the boundary conditions, we have 


u(x, y, t) = 0, v(x, y, 4) =0 on B 


vs 
on 


v(x, y» t) - 0, = 0 


The vorticity is given by? 
f(x, y, #) = — 2dy. (S) 
According to the Navier-Stokes equations the vorticity must satisfy the equation 
og ot 


u— +v0— = vA (6) 
ot Ox 





* Received March 1, 1947. This paper was presented on February 24, 1947 at a joint Colloquium of 
the Department of Physics and the Graduate Division of Applied Mathematics of Brown University. 

1 What is meant here by the term regular curve is only that B fulfills the conditions required for the 
use of the Green formula. 

2 Ay indicating the Laplacian: Ay = 0% /dx2+0*y/dy*. 
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N 


in which v is a constant, the kinematic viscosity of the fluid. Substituting the values 
(1), (2), and (5) into (6), we obtain a partial (non-linear) differential equation for 
¥(x, y, t), which is characteristic for the two dimensional flow of an incompressible 
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viscous fluid. We shall only deal with regular flows, that is with stream functions such 


that the functions 
0¢ o¢ Op ou Ov 0¢ 
Pati => i es 4 = : 
Ox oy Ot ot ot ot 


o- ’ —-s a 


are continuous in x, y, t for (x, y) in D+B and 4, StSh. 

This paper is divided in three parts. In the first part, the Fourier transforms 
Ww, we, t), U (a, we, t), V(w1, we, t), Z(w1, we, t) of the stream function y, the velocity 
components 4%, v, and the vorticity fare introduced and it is proved that ¥, Uand V’ 
are obtained from Z by multiplication by simple rational functions of w,; and w2. Next 


the kinetic energy of the flow 
. 1 \ ~ 
E = ° f (u? + v*)do, (do = dxdy) (7) 
& D 
and its spectral decomposition 
E= f 7(w1, we, t)dw, (dw = dw dwe) (8) 


are considered and it is shown that the spectral function y(w1, we, t) has also a very 
simple expression in Z. In other words, the entire harmonic analysis of the flow can 
be based on the Fourier transform of the vorticity. 

This part of the paper uses only the definitions of y, u, v, § and the boundary 
condition; the results are valid whether the fluid is viscous or not. 

In the second part Eq. (6) is used and some inequalities are given which must 
be satisfied by the Fourier transforms when the variation of the time ¢ is consider 
Some of these inequalities concern the variation of the spectral function with time. 
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All inequalities are based on the fact that the kinetic energy of the flow is decreasing more 
rapidly than an exponential function of t. 

Next, an integrodifferential equation for the Fourier transform Z of the vorticity 
is derived from (6). This equation might be the starting point for the study of the 
variation of the spectral function with the time ¢. 

Finally, in the third part, the case that the flow fills the entire plane 

M4 —-axa<crxr<c4+o, —-~o<cy<+ea 


is considered as the limit of the flow in a bounded domain, the fluid being at rest at 
infinity. The integrodifferential equation simplifies considerably in this case. 
If f(x, y) is a real function of the real variables x and y defined* and continuous 
in D+B, its Fourier transform is given by 
" 1 7 1 
F(w, we) = f(x, yen thot w2wdg, (9) 
Tr D 
where the frequencies w, #2 are real. We shall use the following well-known properties 


of the Fourier transforms.‘ 
a) F(w@;, #:) is a complex function of the two real variables w;, w2 defined in the 


entire plane 
—andw< +o, — 2 Cw< + a, 
b) F( —@1, — 2) = F(w, We), 


‘~*~ 
~ 


where F denotes the conjugate of F. 
c) F(w,, @2) is a continuous function of the two variables a, w2 in every point of 2. 
d) F(a, we) is bounded in Q: 


1 
| F(w,, we) | < if f(x, y) de. (10) 
7 D 


We shall also use the bound 


\ Ss 1/2 
| F(w:, w2)| S - 1 f f(x, y)de| , (11) 
4r? D 


where S=/pdo. Equation (11) is deduced from (10) by means applying Schwarz’ 
inequality. 

e) If | co, | + | we| —+o, then F(a, w)—0 (Riemann’s theorem), 

f) If d(x, y) =af(x, y)+bg(x, y), where a, b are constants, and f(x, y), g(x, ¥) 
continuous functions in D+B having the Fourier transforms F(a, we), G(@1, @2), 
then d(x, y) has the Fourier transform ®(a, w2) =aF(@1, @2) +bG(w1, we). 

If f(x, y) and g(x, y) are continuous functions in D+B having the Fourier 
transforms F(w;, w2) and G(@;, w:), respectively, and if 


F(a, We) = G(w1, ws), 


n order to apply directly most of the known results of the theory of the Fourier transforms, it is 
ften useful to consider f(x, y) as defined in the entire plane X with f(x, y) =0 at every point outside 


+B, 


*S. Bochner, Vorlesungen itber Fouriersche Integrale, Leipzig, 1932, pp. 183-197. 
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then 
f(x, y) = g(x, y) 
in every point of D+B. 


h) If F(a, we) is absolutely integrable, 


[ | F(w1, we) | dw < + , 


the integral 
ke w)etoizto2y) day 
defines a continuous function of x, y for all values of x, y;in D+B: 


f(x, y) = f F(w1, weet @1*+#2V) dy, 


while 


0 = f F(a, woe @1*+¢2 dy 
2 


outside D+B. Since the integral is continuous for all values of x, y, the function 
f(x, y) must vanish on B. Thus, if the given function f(x, y) does not vanish on B, 
its Fourier transform cannot be absolutely integrable.5 

z) If f(x, y).has continuous derivatives 0f/dx, 0f/dy in D+B (which is always the 
case for the functions considered here), one has 


f(x, y) = f Fe. we)et@irto2u) dy (12) 
in every point of D, the integral being now an improper integral defined as the limit 
lim f F(a, weet 12 to2V) dy, (13) 

ate Jo, 


C, being the circle wi+w Sd? (Cauchy’s principal value). As a rule, (12) does not 
hold on the boundary B. 

j) f(x, y) and g(x, y) being continuous functions with continuous derivatives in 
D+B and F(a, w) and G(w;, w.) their Fourier transforms, we have 


i f¢ _ 
- | f(x, y)g(x, yje~*2402dg = f F(w1, w2)G(wi + 41, we + 82)du, (14) 


4° D 


the meaning of the integral being the same as in (12). Here, 6, and 62 are arbitrary 
real variables; in particular for 6; =62,=0 one has Parseval’s formula: 


5 For instance, taking for D+B the square -—1S*x5+1, —1SyS-+1 and assuming that f(x, y) =1 
in D+B we obtain the Fourier transform F(w;, we) =(1/7%wwe2) sin w sin we which is not absolutely 
integrable. Eq. (12) holds, however, with the definition (13). 








ne Sree ere 
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1 = 
f f(x, y)g(x, y)do = f F(w1, w2)G(wi, w2)dw. (15) 
4n’ J p 2 


1. Let us introduce the Fourier transforms of the stream function, the velocity 


components and the vorticity: 


1 . 

V(w1, we, t) = ri =f v(x, y, te *irtewdg, (16) 
tT! D 

1 ; 

U(w1, w2, f) = : f u(x, y, De iiztowwdg, (17) 
mr’ J p 
1 

V(w1, we, t) = r f v(x, y, de eiztowdg, (18) 
7!’ D 
1 

Z(w1, we, t) = ; f t(x, y, Dem *orztew dg, (19) 
TJ p 


The first part of the present paper deals with the purely kinematical significance 
of , u, v, ¢, all computations being supposed to be made at a given time ¢. Thus, there 
is no particular need for stressing the particular value of t. For the sake of brevity we 
shall therefore write W(x, y), V(w1, we), etc. for P(x, y, t), V(wi, we, t), etc. 

THEOREM I—The Fourier transforms of the stream function and of the velocity 
components are expressed in terms of the Fourier transform of the vorticity by means of 





2 

V (a1, we) =_— — Zw, we), : (20) 
wt + wf 

a 2iwe © 

U(w1, w2) = ——— Z(w1, 2), (21) 
wi + we 

; —_ 2tw} 

i (w, w») =- z —Z(w, We). (22) 
wt + wi 


Substituting (1) into (17) we obtain 


1 d 
U =- f =. i(wizto2dg 
4r’ D oy 


1 0 : 1W2 ’ 
nos [pe-iorzteo) |g + - — ve i(wirtoawdg, 
F 4nr? J p 


4r? ) oy 


Application of Green’s formula to the first integral yields 
1 
U=- Bye-iurztowds + jwsW. 
4a? B 


On account of the boundary condition (4) the first integral equal zero. Thus 
U = iwoV. (23) 


Similarly 


V = — ia. (24) 
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Next, let us apply Green’s formula 


ad Of 7] 
f (fAg — gAf)do = — f ( — ‘) ds 
D B On on 


setting f=y, g=e—*@*++, On account of the boundary condition (4) we have: 


S 


f [pAc~ierztorv) — Ape iwirt+ow |dg = 0, 
D 
If we replace AY by —2¢ in accordance with (5); and remark that 


2 2 
Ae i(wiztwey) — — (w; +: wo)e i(wiztwey) (25) 


we obtain 
— (wi + we) Ye ieizto2wdg +- 2f Ce iteiztowdg = () 
D D 


which is Eq. (20). Inserting the value (20) of V into (23) and (24) we immediately 
obtain the expression (21) and (22) and our theorem is proved.® 
Combining (21) and (22) we obtain the following relation 


21Z = w2U — ayV. (26) 


THEOREM II—The spectral function y(w, we) of the kinetic energy E which 1s de- 
fined by (8), is. expressed in terms of the Fourier transform of the vorticity by means of 
; | Z(w1, we) |? 7 

y(w1, we) = 81? - om —— * (27) 


wi + wi 


To prove this, we start from Green’s formula 


Oy” oy? oy 
f pAyda +f (; )+(; ) |e = -f y — ds. 
D D Ox oy B on 


The integral over B equals zero on account of the boundary condition (4). If we re- 
place Ay by its value (5) and dy/dx, dy/dy by (1) and (2), we obtain 


oe f yide. (28) 


Applying Parseval’s formula (15), we find 


E = 49° f WZdes, (29) 


and hence 
¥(w1, @2) = 49° (ws, w2)Z(w1, we). (30) 


6 It is perhaps worthwhile to stress that the boundary condition (4) plays an essential part in this 
proof. As a rule, the relation /j=i#;F does not hold between the Fourier transforms F(w, w2) and 
F,(w1, w2) of f(x, y) and df/dx. For instance for the function f=1 considered in remark 5, F,=0, while 


9 


to F = (ir? /wo) sin w, SiN wo. 





ME Cat ae 
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From this we deduce (27) by inserting the value (20) for WV. 

The theorems I and II reduce the harmonic analysis of the flow to the study of the 
Fourier transform Z of the vorticity. We shall now establish some properties of this 
function. According to (12). 


t(a, y) = f Z(w1, we)etoirtorn) dy 
2 


= lim f 2%. wet @irtway) dyy, (31) 
Cr 


A+ 


We can then interpret Z(w, w2)dw;da, as the contribution of vortices, the frequencies 
of which are between @; and @,+da;, and we and we+da». 
Z(w, #2) is continuous and bounded in the whole plane Q. Noting that 


1 
v(0, 0) = —f vdeo 
4r? J p 


is finite, we obtain from (20), 


Z(0, 0) = 0. (32) 
Equation (26) then yields 
2|Z| s|we||U| + | ol) Vv]. (33) 
On account of (11) we have 
i\U|s : =| f wie] | Vis =| f vie] 
49? D 4r? D 
and consequently 
re S _ 
oe eT eo (34) 


Equation (33) gives then the following upper bound for | Z| : 


2S _, _ 
iS (lor!) +) wo! VE. (35) 
82" 
At this point it seems interesting to note a property of the Fourier transform Z. If 
o(x, y) is an arbitrary harmonic function in D, and ®(a, we) is its Fourier transform, Z 


is orthogonal to ®, i.e., 
ff 2s 2) ®(w, w)dw =s (0, 
2 


In fact, the vorticity ¢ is orthogonal to any harmonic function ¢ in D.’ To prove this 
we need only use Green’s formula 


7 J. Kampé de Fériet, On a property of the Laplacian of a function in a two dimensional bounded domain, 
t derivatives of the function vanish at the boundary, Mathematics Magazine No. 2, 1947. 





when the first a 
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Oy 0g 
f (pAy — pAd)do = — f (+ — ~) ds. 
D B On On; 


On account of the boundary condition (4) the second integral vanishes. Setting 


Ay = — 2, Ag = 0 


f otdo = 0 
D 


which, on account of (15), is fully equivalent to our statement. 
Let us now consider the spectral function y(a1, w2). We can interpret 7(w1,We)dw,dwe 
as the amount of kinetic energy coming from the vortices with frequencies between 


in the first integral we find 


o, and w, +d, and w. and we+dan». 

The Fourier transforms V and Z being continuous and bounded in the entire plane 
2, we see from (30) that (w:, w2) is also continuous and bounded in the entire plane. 
Moreover, y(@1, #2) 20 for all values of w, and a; on account of the same formula 


¥(0, 0) = 4x20, 0)Z(0, 0), 


and hence 
(0, 0) = 0. (36) 


More precisely, if we write Eq. (30) in the following manner 


(w1, w2) = 2n?(ws + we) | W(wr, ws) [2, 





we see that 


¥(w1, we) i 


v(0, 0) |2. (37) 





lim 2x? 
@,0,4,-0 wt + w} 
Since V and Z tend towards zero when |«:| +|«:|->+ ©, we have 


lim y(w1, we) = 0. (38) 
lw, |+]w.|—--+ 2 


From (27) it follows then that y(w:, w2) is decreasing more rapidly than 1/(wj+w). 
Combining (27) and the inequality (35) we obtain an upper bound for y(a1, 2): 


S (lor| +] eal)? 





y(w1, we) S > ’ 
4r? wi + we 
and hence 
_ 
y(w1, 2) < ae E. (39) 
23? 


2. Let us now consider the variation of the flow with the time ¢, in accordance 


with Eq. (6) for a viscous incompressible fluid. 

We shall use the following proposition :* The kinetic energy is always decreasing 
8 J. Kampé de Fériet, Sur la décroissance de l'énergie cinétique d'un fluide visqueux incompressible 
occupant un domaine plan borné, C. R. Acad. Sci. Paris 223, 1096-1098 (1946). 
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(dE/dt <0); more precisely, it is decreasing more rapidly than an exponential function of 


Sv 
E = Eo exp ( — ‘). (40) 
V\ K D 


Here, Ey is the value of E at the time t=0, and Kp is a positive constant depending 


, 4 


the time: 


only on the domain D. 

Applying this result to (35) and (39), we have the following 

THEOREM III—The upper bounds of | Z| and ¥ given by the inequalities (35) and 
(39) are always decreasing, more precisely, 


, \ 2SEo | 4p 
Z(w1, we, t)| S$ — (| w:| +] w2|) exp( — —=-2 (41) 
| 8x? VKp 
SEo ( Sv ) (42) 
"y(wW), Wo, t) = are exp —_ ae — pe | 
2x? ;' V/Kp 


A more concrete picture is furnished by the surface~y(w, we, 4). This surface passes 
through the origin O and is asymptotic to the plane y =0 for large values of a, we; 
its possible maxima are all below a given plane y = MV; with increasing time this plane 
approaches the plane y =0 in an exponential manner. 

In order to study the variations of the spectral function (a1, we, +) with the time 
i, we shall establish the equivalent of Eq. (6) in terms of Fourier transforms. To this 
end, we must compute the Fourier transforms of 0¢/dt, of ud¢/dx+vd¢/dy and of Af. 





Fic. 2 
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For the first function the result is obvious. Since ¢ and 0¢/dt are continuous in ¢, 
the functions Z and 0Z/d¢t are also continuous in ¢. Thus, 


Go. 1 0¢ 
— Z(w}, w2, t) = - f — etlirtewdg, (43) 
ot 4r? J p A 


The Fourier transform of 0¢/dt is therefore given by 0Z/dt. 


To compute the second Fourier transform, we start from the equation 


o¢ 0¢ 2: re) 0 
u - + y - e7 ttwiztesy) = — [wce-i(orrtoay) | a mae [vge-i(orz+o2w) | 
Ox 


Ox OV oy 
+ tail + woevje oirtory), 


We integrate this equation over D and note that by Green’s formula, 


0 
f — [uge-irzte |dg = -f auveMoiztervds = Q), (44) 
“ D B 


Ox 


0 
f _—_ [ote- i wiz twa) |dg = — f Bute—ierto2u)ds ——— 0, (45) 
D oy B 


on account of the boundary condition (3). Accordingly, 


o¢ 0¢ ; 
“u— of vy — pe ~tiztonwdg = iw) ucem' wittwov) dg 
D OV Oy D 


+ iws f veetoiztoon) dg, 
D 


Using formula (14), we obtain 


1 a 
- J uvetwirte2n) dg = f U (41, 62)Z (8; =} @1, A. oy w2)d8, 
7” D P 6 
1 ; ns 
2% f vveMoirtewdg = f V (02, 02)Z(01 + wr, 02 + we) dd, 
4r*J p =) 


where the variables of integration are now called 6,, 62, and d@ denotes d6,d@2, the inte- 
gral being extended over the entire plane 0. We thus have proved that the Fourier 
transform of ud{/dx+v0¢/dy equals , 
if [woiU (61, 62) + weV (04, 62) |Z (01 + ®1, A. + we)d6 
e 


or, what is the same by (21) and (22), 


v 
Awe — Bow s 
2 Te 2 Z(1, 62)Z (4 + w1, 62 + w2)d0. (46) 
e 6+ 43 
Let us now consider the Fourier transform of the Laplacian of the vorticity: 


1 
C(w1, we, 1) = if Ate i@izteavdg, (47) 
D 


4? 
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Applying Green’s formula 


J [Ager tesetons) — gale Merste20)) Jo 
D 


Ta a 
= -f — e~ i(wiztwry) — t— (e~ *(e1z+e2y)) ds 
pLOn on 


and setting 


1 0 
(w, we, 1) = —f [= + i(wra + o:8)r| e~ Kmotuends, (48) 
4n? J 2 Lon 
we obtain 
C(w, We, t) = (0; + ws)Z (wr, We, t) aad d P(w, We, t). (49) 


Reviewing the expressions of the Fourier transforms of 0£/dt, u0{/0x+-v0f/dy, and 


At we may state the following 
THEOREM I \V—Equation (6) is equivalent to the following integrodifferential equa- 
tion for the Fourier transforms: 





0 4 Aw. — Bow = 
- Z(w, We, t) aad 2f aa a -Z(6;, 60, )Z(0; + W1, 6 + Woe, t)dé 
ot . 6; + 4 


=— »(c0s ae wnZ (en, we, 1) — vb(w), we, 2). (50) 


This equation, being fully equivalent to the Navier-Stokes equations, is the rational 
starting point for any rigorous study of the Fourier transform of the vorticity and 
hence, of the spectral function. This study seems to be very difficult; perhaps it will 
be only possible to check, in some approximate way, some of the working hypotheses 
which were assumed in recent papers, e.g., the assumption that the big eddies have a 
tendency to degenerate into smaller ones. This means that when ¢ is increasing the 
peak on the spectral surface y(wi, w2, t) is gradually shifted towards the high fre- 


QueNCies @, We. 
3. We can get interesting general results, if we consider the case where the flow 
extends over the entire plane. 
X —-a <c<x<+o, —a<c y+ o. 


This case may be regarded as a limiting case of previous problem. The boundary 
condition (3) means now that the fluid is at rest at infinity: 


lim u(x, y, 4) =0 lim v(x, y, 4) = 0. (51) 
r=+tf/x?+ y2. 


If the kinetic energy is to remain finite, 


1 
E= - f (u? + v)da < + ~, (52) 
< xX 
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the decrease of u and v must be sufficiently rapid. A sufficient condition is 
u = O(r-3!*), v = O(r-*/*) for large values of r. (53) 


To insure the existence of the Fourier transform 


1 
Z(w1, w2, 4) = - - foc y, De t@iztowdg, (54) 
4n? J x 
it is sufficient to assume that ¢ is absolutely integrable 
fi ltlac< +o. (55) 
x 


In this case, all the properties a) to j) enumerated for the Fourier transform (9) still 
hold for (54). The condition is obviously satisfied if: 
¢ = O(r-*) _ for large values of r. (56) 
To compute the Fourier transforms of 0¢/0t, ud¢/dx+vd¢/dy, and Af, we can use the 
same process as before, taking for D the circle x?+y?<R? and then letting R tend 
towards -+- ©. We must observe, however, that we no longer have u=0, v=0 on the 
boundary B; we must therefore carefully examine some terms in our equations. 
In the evaluation of the Fourier transform of ud{/dx+v0¢/dy, the two terms 
(44) and (45) do no longer vanish. According to Schwarz’ inequality we have, how- 


1/2 1/2 
ae: ee en ri .«t . 9 
f auce twirtwsy as | <{ | u||¢jdss lf v's | f rs | . 
B B B L774 


From our assumption concerning the decrease of u and §, it is obvious that 


lim f u’ds = 0, lim f ev*%ds = 0. 
Rate JB R++o JB 


Similar considerations apply to (45). Thus, at the limit, the expression (46) for the 
Fourier transform still holds. 
In the evaluation of the Fourier transform of Af we obtain a very fortunate simpli- 


fication: from (48) ,we have 


fgihcias 1m. <5 Pr ee ee ee 
| ® |=) » le + ( | @1| + | @2| ) | $ I ds. 
4r? JL) On 


The assumption of decreasing of { yields 


ever, 


O¢ | : 
lim f | ——|ds = 0, lim | ¢| ds = 0; 
R-+= J 3|dn| R-+o JB 
thus 
lim @=0. 
Ro+« 


The Fourier transforms C and Z being extended to the whole plane, we have thus 


now : 
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CcC=-—- (co, + e0s)Z. 


THEOREM V—In the case of a flow extending over the entire plane X , where the fluid 
is at rest at infinity, and where u, v, € decrease so rapidly that 





1 
e-—f (u? + v*)\da << + 
- 


and 


filtla<+e, 
x 


the Fourier transform of the vorticity satisfies the integrodifferential equation: 


rs : O1wW eS 60w om 
Z(w, we, t) + 2f —.—, — Z(81, 92, t)-Z (01 + w1, 02 + we, é)d0 (57) 
at e 6;+ 6 
2 2 
= — y(w: + w2)Z(w1, we, Z). 
This equation seems to me the most simple that can be obtained in this manner. 


To give at least an application, we shall establish the general solution of (57) for 
the particular, but interesting, class of flows, where the vorticity is constant along 
every stream line. 

The flows characterized by {=F(W) have been the subject of interesting re- 
searches.® It is in this case and only in this case that non-linear terms disappear in the 
Navier-Stokes equations; their Fourier transform vanishes, of course. Equation (57) 
is then reduced to 

az 2 2 
— = — ro(w + w,)Z, 
ot 


and the general solution is 
2 2 

Z(w1, @e, t) oa Zo(wr, woe” orton) t, 
Here, Zo(w1, w2) denotes the arbitrary value of Z for ‘=0. It is seen that, in this case, 
not only an upper bound of Z, but Z itself is decreasing exponentially for every value 
of Wi, We. 
® Ratib Berker, Sur quelques cas d’intégration des équations du mouvement d'un fluide visqueux incom- 
pressible, Thesis, Lille, 1936. 
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THE METHOD OF THE HYPERCIRCLE IN ELASTICITY 
WHEN BODY FORCES ARE PRESENT* 
By 
J. L. SYNGE 
Carnegie Institute of Technology 
In an earlier paper! the method of the hypercircle in function space for boundary 
value problems in elasticity was developed under the restriction that body forces 
are absent. The purpose of the present note is to remove that restriction, and present 


the method in complete form, using the notation and basic ideas of PS. 
Let an elastic body (in general anisotropic) be subjected to body forces X; per 


unit volume, so that the equations of equilibrium are 

Ei;,4 + Xi = 0. (1) 
As regards boundary conditions, we shall suppose that the surface B of the body may 
be divided into regions designated by 


[in, te], [utn, T |, [, Tn], 7. T |. (2) 


Here [u,, u:| indicates a portion of B on which both the normal component 4, of 
displacement and the tangential (vector) component u,; of displacement are assigned ; 
[u,, Z:| indicates a portion of B on which the normal component of displacement 
and the tangential component of stress are assigned; and so on. 

Asin PS, a state of stress throughout the body corresponds to a vector in function 
space, and the scalar product is defined by S’-S’’ = fe/,E{/dv, the integral being taken 
throughout the body. We set up the following notation: 

S=natural state satisfying the equations of equilibrium and compatibility, 
and also the boundary conditions. 
S’ =basic stress state, satisfying the equations of equilibrium 
Ei;,;+ Xi = 0, (3) 


and all boundary conditions on stress; this means that 


Ti =T, on [tp, T |, 
Ti =T, on [ we, Fal, (4) 
Ti =T, and Ti =T, on [T,, T:]. 

I/ (p=1, 2, - - » m)=orthonormal homogeneous stress states, each satisfying the 


I 


equations of equilibrium without body forces, 
_/ ~ 
E (iii = 9, (S) 
and making 7%,), vanish wherever T, is assigned and making T(,, vanish 
wherever 7; is assigned. 
S’’ =basic displacement state; this corresponds to a state of stress which satisfies 
* Received Aug. 25, 1947, 
1W. Prager and J. L. Synge, Quarterly of Applied Mathematics, 5, 241-271 (1947); this paper will 
be referred to as PS. Cf. also J. L. Synge, The method of the hypercircle in function space for boundary value 
problems, Proc. Roy. Soc. London (A) 191, 447-467 (1947). 
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the equations of compatibility, so that a displacement u/’ exists. This dis- 
placement satisfies all boundary conditions on displacement, so that 


u<’ =U, and «;’ = 4, on ttn, u,|, 
u,’ = 4%, on [ 2¢,., T. |, (6) 
ui’ =u, On [ wr, mt 


responding to a state of stress satisfying the equations of compatibility, so 
that a displacement (J, exists; this displacement is such that u“/(/), vanishes 
wherever u, is assigned and u(j, vanishes wherever u; is assigned. 

It is important to note that, if k’ and k’’ are any constants, then S’+’I; be- 
longs to the class of basic stress states (S’) and S’’+k’’I?’ belongs to the class of 
basic displacement states (S’’). 

We now proceed to prove that the extremity of S lies on a hypersphere determined 


by S’ and S’’. By Green's theorem we have 


- oa iMunt in 
Ss’-S = f eisBido = fui, Biss 
si 7 
= f uy reap — fut Bis to, 


wnere dB is an element of the bounding surface B. Hence, by (3), 


S’.S” = fuer aB+ ff ui*Xede. (7) 


I/’(q=1, 2,-++mn)=orthonormal homogeneous displacement states, each cor- 


But S belongs to the class of S’ and also to the class of S’’. Thus we can substitute 
S for S’ or S”’ in (7) and so obtain the following equations: 


fue riap + f ux, 


S’-S = f arias + f wXede, (8) 


S-S = f wrap t f uXede 


By addition and subtraction of (7) and (8), we obtain 


S:S” 


II 


S:S-S-S’ —S-8”"+S'-S” =f — ui’)(T; — Ti)dB. (9) 


For any pair of vectors, U;, V;, we have 


f vwas = f v.v.dB +f ura, (10) 


where U,, V, are normal components, and U,V; is the ordinary scalar product of the 
tangential vector components. We apply this reduction to the right hand side of (9), 


ame 
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and then split each integral into four parts corresponding to the four regions (2). It is 
easily seen from (4) and (6) that the result is zero, and so (9) reads 
S-:S —S-S’ —S-8”+S’-S” = 0, (11) 
or, equivalently, 
[Ss — 3(8’ + 8”)? = i(S’ — 8”)*. (12) 
This formula locates the extremity of S on a hypersphere in function space with 
center at 
Cy = 3(S’ + 8”) (13) 
and with radius R given by 
R* = 3(S’ — 8”)’. (14) 
Since all that is required of S’ and S”’ is that they belong respectively to the 


classes of basic stress states and basic displacement states, it follows from a remark 
made earlier that we may make the following substitutions in (11): 


for S’ substitute S’ + &’1/, 
for S” substitute S” + k’Iy’, 
where k’, k”’ are arbitrary constants. When we do this, and then subtract (11) from 


the result, we get an equation which must be satisfied for all values of k’, k”’. This 
vields 
S-I, = §”.I/ (p= ‘ Whe: (15) 
S-ly’ =S'-ly’ (q=1,2,---,m), 
and also 
I, -If’ = 0 (p= 1,2,---,m;q=1,2,---,m); (16) 

this last result shows that the set of homogeneous stress vectors is orthogonal to the 
set of homogeneous displacement vectors. 

Equations (11) and (15) locate the extremity of the natural vector S on a hyper- 
circle T’. It is easy to prove that the center C of I’ and its radius R are given by 


1 m n 
7 Is +8”"- DUE {s' -s")-L} + UU {s' - ”)-13"} |, (17) 
é p=1 q=1 
1 m n 
R? = r | s —s”)?-— > {(s’-—8”)-1/}?- > {S’ -— 8”)-Ly’ +]. (18) 
p=1 q=l 
Let us use I, (r=1, 2,--+,m-+m) to denote the whole set of homogeneous 


vectors I/, I/”, so that, by (16), I, form a set of orthonormal vectors in function 
space. The equation of the hypercircle ! may be written 


S=C+ RJ, (19) 
where J is arbitrary except for the equations 


Ji=1, . J L=0 (r= 1,2,---,m+n). (20) 
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We now seek the points on I at minimum and maximum distances from the 
origin O of function space. For any point S on T we have, by (19), 


S? = C?+ R?+ 2RC-J, (21) 


and so, on account of (20), for stationary values of S? we have 
mrn 
C=aJ+ > dl, (22) 
r=} 


where a and b, are undetermined multipliers. We deduce 
b, = C-L,, (23) 


and 


m+n 


a? = | es 1(C-1) | = c:— > (C-1,)% (24) 
r=1 r=1 
For a stationary value of S? we must take, by (22) and (23), 
J=a Le -> 1(c-1) |. (25) 
r=] 


We find it convenient to define a vector G by 





G=  s/ +S” - = I,{(S’ + s")-13], (26) 
£L r=1 
so that 
oe = “ is +8”)? — 3 {(S’ + s")-1}3|. (27) 
By (17) and (26), we have 
C=G+ > IZ(S”-1f) + i 1,’ (S’- Il’). (28) 
p=1 q=1 
We note that G-I,=0, and so by (28) 
C-G=G’*2 0. (29) 
By direct algebra we may prove that 
C~ F464) «6. (30) 
r=1 


Thus, by (24), a= +|G], and, by (25), J=+ G/| G|. By (19) the stationary values 
of S?, if S ranges over I’, correspond to the vectors 

C + RG/|G|. (31) 
Of these two vectors, that with the + sign is the greater, on account of (29). Accord- 


ingly, we may say that the point of I farthest from the origin of function space has 


the position vector 
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V’ = C+ RG/|G}, (32) 
and the point nearest to the origin has the position vector 
Vv” =C — RG/|G|. (33) 
Thus we have the following inequalities for the magnitude of the vector S correspond- 
ing to the natural state of the elastic body: 
v2 ¢ ¥", (34) 
or, by (29), 
C? + R? — 2R|G| < S? < C?+ R*+4 2R|GI, (35) 


where C is given by (17), R by (18), and |G| by (27), on taking the square root. We 
can show by direct algebra that 


1 
2 


| = {(S”-1;)? — (S-1,)'] 


p=1 


1 
C? R? = = (S’? + S’’?) 4 


+2 {S'-l’)?- "44 |. (36) 


q=1 
The inequalities (34) or (35) bound the strain energy of the natural state below and 
above. It is interesting to note that the usual minimal principles bound the potential 
energy of the natural state, rather than its strain energy. 
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GENERAL IMPEDANCE—FUNCTION THEORY* 


BY 
PAUL I. RICHARDS 
Brookhaven Natignal Laboratory 


1. Introduction. In the field of lumped-constant circuits, considerable theory is 
available’ concerning relations satisfied by all physically realizable impedances. Condi- 
tions which are both necessary and sufficient for physical realizability have been de- 
rived. It is the purpose of the present paper to extend some of these techniques to more 
general, distributed-constants circuits. 

The basic method of existing theory consists in extending the impedance Z = R + 
iX to “complex frequencies” s = y + iw so that Z becomes a complex function of a 
complex variable. In the field of lumped-constants, the circuit functions are rational 
and the extension is effected merely by substituting s for iw. 

2. The concept of impedance. When we pass from the field of lumped-constant to 
more general circuits, we can no longer say that the impedances are rational functions 
of w. This circumstance raises a number of difficulties. First, we must consider singu- 
larities other than poles, for it is known that all non-rational functions have at least 
one singularity which is not a pole. (For example, the input impedance of a shorted dine, 
tanh (iwl/c), has an essential singularity at infinity.) Secondly, we no longer have the 
relatively simple type of differential equation which occurs in lumped-constant circuits. 
At our level of generality, there are available only Maxwell’s Equations. These con- 
siderations make it advisable for us to review anew the general concept of impedance 
and the means of extending this concept to complex frequencies. 

Before embarking on our general discussion, we must stress a point frequently over- 
looked in the foundations of electro-magnetic theory. When introducing this subject, 
the remark is often made that specialization to harmonic time dependence is of no great 
consequence due essentially to the Fourier Integral Theorem. At a later stage, it is 
then mentioned that the material constants in the resulting equations may depend 
upon frequency. Often it is not mentioned that, in the case of such dispersive media, 
it is obviously impossible to return to the original Maxwell’s Equations for arbitrary 
time dependence. 

Thus, if we wish to be accurate, we should not say merely that the postulates of 
electro-magnetic theory are contained in Maxwell’s Equations. Rather, we must say, 
at least for dispersive media (and, in actual practice, in general), that the postulates 
consist of: (a) Maxwell’s Equations for harmonic time dependence and (b) Fourier (or 
Laplace) superposition for arbitrary time dependence. 

Our arguments will be based essentially upon (b). The purpose of the discussion 
was to emphasize that (b) is a basic postulate of electro-magnetic theory and is not a 
mere pseudo-mathematical consequence thereof. 


*Received May 7, 1947. The results and general method of this paper conform to a section of the 
author’s doctor’s thesis, the research for which was carried on under the direction of Dr. Ph. LeCorbeiller 
at Harvard University. The arguments, however, have been markedly changed from those originally 





employed. 
1For an excellent summary of lumped-constant impedance-function theory, see H. W. Bode, Network 
analysis and feedback amplifier design, D. Van Nostrand, New York, 1945. 
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We must next decide what we shall mean when we speak of an impedance in general. 
We certainly shall want any impedance Z to satisfy a relation of the type: 


(response quantity) = (applied quantity) K Z(iw). (1) 


Actually, we shall place no restrictions upon the quantities involved except for those 
limitations already implied by the form of (1). These implied restrictions, however, 
merit considerable discussion. 

Firstly, note that applied and response quantities cannot always be interchanged. 
A familiar example in the field ot lumped constants would be generator-voltage and 
short-circuit-output-current for a four terminal network. Under no conditions can the 
latter be considered “applied”. (It is this impossibility of interchange which causes the 
analytic properties of a transfer function to differ markedly from those of its reciprocal.’) 

Secondly, it is specifically stated in (1) that Z is to depend only upon frequency. 
This implies that (z) the circuit configuration is constant in time and (72) the circuit is 
linear (i.e. the material constants depend at most on frequency and spatial position). 
We shall add the further restriction that (ii) the circuit is passive; mathematically, 
this is the statement that the total power absorbed is eventually non-negative. 

Finally, (1) implies that the applied and response quantities can be characterized 
by two single numbers, one of which is determined by the other. Thus the circuit must 
be so constructed that the entire spatial dependence of all field components is com- 
pletely determined by frequency, while the absolute magnitude of the field at all points 
is determined by a single “excitation parameter’’. 

We shall not attempt an exhaustive discussion of the cases in which this fourth 
restriction can be met, but a few remarks are in order. The first picture that comes to 
the reader’s mind when this requirement is stated is probably a circuit which—while 
possibly containing waveguides, antennas, etc.,—has definite small input terminals. 
But what does “small” mean when the frequency has not been specified? Consider 
lumped-constant impedance-function theory. Here, definite input terminals are assumed 
and yet the impedance functions are treated as having physical significance beyond 
the frequency of visible light and, indeed, at w = ©. What is essentially done, of course, 
is to idealize the physical picture in a self-consistent fashion so as to obtain a theory 
which is tractable and, at the same time, physically significant within the frequency 
range of interest. It seems reasonable that the same procedure could be carried out 
over greater frequency ranges; a case in point is transmission-line theory, which avoids 
the lumped-constant approximation of treating a line merely as a pair of leads but which 
ignores the effects of radiation, waveguide modes, and discontinuity reactances. 

Following these thoughts one step further, we note that we may dispense with the 
concept of definite input terminals. We might, for example, choose as the quantities in 
(1), the transverse dominant-mode components of the forwardly-propagated electric 
and magnetic fields at some point of a reference plane located within an input wave- 
guide. The characteristics of the waveguide would assure satisfaction of the above con- 
dition provided that the higher modes could be excited only through the dominant 
mode. This could be accomplished by specifying that the guide be fed by an input 
guide which “‘tapers to a point generator’’. In practice this last phrase could mean that 

*Note that the usual roles of applied and response quantities have been interchanged. This facilitates 


treatment of transfer impedances. 
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the taper is to be cut off and attached to an actual generator after the cross-section 
has become too small to support higher modes; this would involve a set of generators 
with associated coupling tapers of different lengths for each frequency range (no taper 
is needed below the cut-off of the main guide), but the impedance as defined above would, 
nevertheless, be well-defined. 

Having reached this level of generality, we are faced with the question of distin- 
guishing between transfer and driving-point impedances. We shall adopt the following 
definition, the motivation for which is obvious in the light of the wave-guide example. 

Definition: If the product of the quantities in (1) has the same algebraic sign as the 
total power (or energy) absorbed by the circuit, then Z is called a driving-point im- 
pedance. 

Before closing this section, another consideration should be mentioned. Inasmuch 
as our arguments will be based upon physical principles, we must require that the cir- 
cuit can actually be constructed. The significance of this remark will, we hope, become 
clearer as the argument progresses. For the present, we may remark that a transmission 
line or waveguide of infinite length must therefore be excluded; though Maxwell’s 
Equations are capable of treating this fiction, the device cannot actually be exhibited 
in the world of reality. An antenna radiating into empty space (or into a completely 
absorbing box if we wish to avoid the non-linear ionosphere) can be constructed and used 
and is therefore not to be excluded. Secondly, we should note that truly lossless 
circuits are also a fiction. We do not wish to exclude these; they are far too useful. How- 
ever, recognizing their idealized nature, we shall, when convenient, treat them as the 
limiting case of a lossy circuit. This corresponds to the physical facts: losses are present 
but, in some cases, effects due specifically to dissipation are small enough to be neg- 
lected. 

3. Necessary conditions. We now wish to extend our impedance Z(iw) to complex 
frequencies s = y + tw. Maxwell’s Equations certainly have formally valid mathe- 
matical solutions for time dependence of the form e”‘ cos (wt) = Re(e**). These solutions 
can be obtained by mere substitution. With dispersive media, however, there would be 
considerable doubt as to the correct value of the material constants at such complex 
frequencies. Moreover, although the substitution process will usually lead to an analytic 
extension of the original function of a single variable, there exist functions (e.g. jump 
functions) which cannot be so extended. We shall, therefore, adopt a very different 
approach, which allows the use of the powerful theory of the Laplace Transformation. 

Let us cause the applied quantity, f(t) say, to vary with time in an arbitrary manner, 
subject only to the condition f() = 0 for t < 0 (without which we could hardly be said 
to “apply” the quantity). Let g(é) represent the response quantity; since the circuit 
cannot respond before it is stimulated, g(t) = 0 for t.< 0. If, for the moment, s = 7a, 
then, by the fundamental superposition postulate of electro-magnetic theory ((b) of 
Section IT), we have: 


[ e**g(0) at 
Z(s) = > nee (2) 
| e"'f(t) dt 


If we now allow s to take on any complex value, the integrals in (2) become the defini- 
tion of the Laplace Transformations of f and g. Thus, by analytic continuation if neces- 
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sary, Z(s) has been defined for all complex s. The superposition postulate of electro- 
magnetic theory implies, in particular, that Z(s) as defined by (2) takes the same values 
on s = iw whatever the function f(t) may be. We shall see later that Z is analytic in 
y > 0° and, at least, on portions of y = 0. By an identity theorem‘ for analytic functions, 
it then follows that Z(s) is well-defined by (2) despite the arbitrariness of f(t). We shall 
also see later that Z has the usual physical significance in y > 0. 

The definition (2) has the great advantage that we may, in particular, choose har- 
monic time variation of extremely low frequency and thereby utilize established facts. 

In accordance with the last remark of Section II, we shail first consider circuits 
of arbitrarily minute but nevertheless finite loss. Moreover, let f(¢) vary harmonically 
in time at such a frequency s, = iw, that the associated wavelength is far greater than 
the overall (“‘metallic’’) dimensions of the circuit. By this statement, we do not mean 
that circuit radiation is to be neglected in the entire analysis but merely that the above 
wavelength is to be so chosen that radiation at that frequency will be extremely small. 
(For.example, a centimeter antenna may be considered as a capacitance at a frequency 
of one cycle/sec.) Thus, though the circuit may in a sense be infinite, the wavelength 
is to be chosen very large compared to its material, physical dimensions.” 

Under these conditions (f(t) = e**', t > 0, 8 imag.),° the (presumably) established 
theories and facts of low-frequency electricity show that the response quantity will be 
given by 

0 t <0, 
g(t) = (3) 
Ae*’’ + h(t) t>0, 


where h(t) represents the transient response and is bounded by e *’ with small but finite 
a.° (Recall our momentary restriction of finite losses.) 
Carrying out the Laplace Transformations in (2), we have: 
Z(s) = A + (8s — %)H(s), (4) 


at 


where H(s) is the Laplace Transform of h(é). Since h(t) is bounded by e~“‘, it follows’ 
that H(s) is analytic in the region Re(s) = y > —a. Thus we may finally state the fol- 
lowing theorem. 

For a circuit with finite loss, however small, there exists an a > 0 such that Z(s) 1s analytic 
iny > —da. 

Setting s = s, , (4) immediately gives Z(s,) = A, the correct value according to (3). 
If we choose another value for s, , the transient excitation may change somewhat (as 
it must by (4) since Z is well-defined), but its general character cannot be greatly altered 
since transient behavior is essentially independent of frequency. (A physicist accustomed 
to working at high frequencies—and thus with short times—may, however, insist on 


’Throughout this paper the phrase “in y > 0’’, in accordance with the usual convention, does not 
assert anything about the point s = ~. As mentioned earlier, an essential singularity often occurs there. 
‘C. Carathéodory, Conformal representation, Cambridge Univ. Press, 1941 (reprinted), p. 82. 

5Note that it is precisely at this point that the fictitious infinite transmission line violates the restric- 
tion of ‘‘physical reality” laid down in Section 2. 

*For simplicity, we use complex notation. The reader can easily verify explicitly that the same final 
results are obtained for the true time variations. 

7G. Doetsch, Theorie und Anwendung der Laplace Transformation, Dover Publications, New York, 
1943, pp. 43, 194, 401. 
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preserving a “‘fine-structure’’ which the audio worker chooses to neglect. For example, 
compare the response to a unit step of a short lossless line terminated by a resistor with 
that computed by assuming the line to be a small inductance—or even a mere pair of 
leads.) Thus H(s) will remain analytic in y > 0 and therefore bounded at any point 
in this region. Consequently Z(s,) will reduce to the correct value A whatever complex 
value of s) is assumed. The usual physical significance may therefore be attached to 
Z(s) for complex values of s lying in y > 0. 

We also note that if s is real in (2) then Z(s) is real and thus, by the Schwarz re- 
flection principle, 


Z(s*) = Z*(s). , (5) 


This also follows formally from the reality of the general form of Maxwell’s Equations. 

We must, finally, consider circuits which, under appropriate conditions and perhaps 
only within certain frequency ranges, may be considered “lossless”. Inasmuch as such 
circuits, for reasons given. éarlier, should really be considered as the limiting case of a 
lossy circuit, we may suspect that the theorem italicized above will “essentially” hold 
in y > 0, and this will prove to be true. 

At sufficiently high frequencies, a dissipationless circuit must be a lossless cavity. 
For lower frequencies we may arrive at the same situation in the following way. It 
may be that the circuit is not truly a cavity but produces practically zero radiation. 
Since this very small radiation is neglected anyway, we may surround the circuit by 
a sufficiently large (but finite) perfectly conducting box; the reflected radiation will 
then be of the sam®@ order of magnitude as that originally neglected. 

Now we wish to re-examine the argument leading to (4) and the theorem following 
therefrom. Certainly we may still choose an s, whose wavelength is very much larger 
than the circuit and its surrounding box. The transient response h(t) in (3) is, again by 
the established facts of low-frequency systems, bounded but never dies away. The 
transient will consist of a superposition of the resonant frequencies of the circuit. Now 
it is shown in the theory of eigen-value problems that, for a lossless cavity, these - 
resonances occur at discrete, isolated frequencies (i.e. they form a sequence with no 


finite limit point). Thus we may say* 
h(t) = > a,e**"*. (6) 


(An added exponentially decaying part may occur if the circuit is not “lossless” at all 
frequencies; this, it is easily seen, contributes no new information about the region 


y > 0.) The Laplace Transform of (6) is 
H(s) = bx a,/(s — tw,). (7) 


Substituting (7) into (4), we see that Z(s) is analytic in y > 0 (as follows merely from 
the boundedness of h(#)),’ and, since the w, have no finite limit point, we may summarize 
our results as follows. These are valid for both transfer and driving-point impedances. 

(a) Z(s) is analytic in y > 0 except only for isolated poles on y = 0; 

(b) Z(s*) = Z2*(s). 

Finally, we wish to find an added condition to characterize driving-point impedances 


8R. Courant and D. Hilbert, Methoden der mathematischen Physik, Interscience Publishers, New York, 
1943, vol. 1, pp. 358, 384. 
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as defined in Sec. 2. This has been done by Brune’ who has shown that if the integrated 
product f(#)g(é) is always positive for all sufficiently large time intervals then: 


Re(Z) = R>0 in the region y¥>0 (8) 


This argument establishes (8) for lossy circuits, and treatment of a “lossless” circuit 
as a limit shows that (8) is valid generally. 

4. Consequences of necessary conditions. Throughout the remainder of this paper, 
only driving-point impedances will be considered. For brevity, we shall refer to such 
merely as “impedances”. 

Any physical impedance must satisfy any theorems derived from the conditions 
italicized above and Eq. (8). The author has made a study of the class of functions 
defined by these conditions; this investigation is purely mathematical and has appeared 
elsewhere.” The results of physical interest will be summarized in this section. 

The only impedance which is purely resistive at all (real) frequencies is the constant 
resistance Z = R, .” 

The purely reactive functions have many properties analogous to those given for 
lumped-constant reactances in the well-known Foster’s Reactance Theorem. General 
reactances are analytic except on the real frequency axis, and there the zeros and poles 
alternate, while the origin is always a zero or pole. The derivative is always positive 
(never zero), and for all w, 


dX/dw > | X/w (9) 


Previously, (9) had been obtained only from Lagrange’s equations for a lumped-constant 
network; it is now known to hold for all generalized reactances. 

A general reactance is necessarily a complicated function, for, if it has only a finite 
number of poles on the real frequency axis, it must be a rational function and therefore 
a lumped-constant reactance. Nevertheless, two expansions used for lumped-constant 
reactances can be extended to general reactances; the processes involved naturally 
become infinite. All of the expansions quoted below converge absolutely and uniformly 
in any bounded region. (Except, of course, for terms involving poles within the region.) 

One such expression is the ‘‘partial fraction” expansion: 


Z = As + (K./s) + >> 2K,s/(s° + w2) (10) 


where w, are the anti-resonant (angular) frequencies and K, , the residues at these 
poles (K, may be zero); A = lim Z/s as s approaches © in| arg s| < 7/2. 
Schelkunoff has suggested’ that there might exist such expansions and, pending 
their proof, has shown in considerable detail how they may be used to calculate the 
input reactance of complicated distributed-constant structures. The anti-resonant 


°O. Brune, Synthesis of a finite two-terminal network whose driving-point impedance is a prescribed 
function of frequency, J. Math. Phys. 10, 191-236 (1931). 

10P, Richards, A special class of functions with positive real part in a half-plane, Duke Math. Journal, 
14, 777-786 (1947). 

Proved for lumped-constant circuits by O. Zobel, Distortion correction in electrical circuits with 
constant resistance recurrent networks, Bell Syst. Tech. J. 7, 438-534 (1928). 

25, A. Schelkunoff, Representation of impedance functions in terms of resonant frequencies, Proc. 
I. R. E., vol. 32, pp. 83-90 (Feb. 1944). 
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frequencies are, of course, relatively easily evaluated. The K,, (i.e. equivalent L and C) 
can be found from energy storage. The first term is less easily calculated; it can be shown 
to be (d/ds)[Z — (K,/s)] at s = O minus the series p ley 2K,,/w. . For this reason and 
also because the series (10) converges rather slowly, it is preferable to employ the fol- 
lowing. Use Z or Y = 1/Z, whichever is zero at zero frequency. Then, setting Z’(0)— 
or Y’(0)—equal to L (the low-frequency inductance—or capacitance), we have 


Z = Ls — 28° & K,/|w(s’ + w3)] (11) 


This series converges much faster than (10) and the terms are more easily evaluated. 
For examples of such calculations, the reader is referred to Schelkunoff, op. cit. 

The validity of (10) justifies the common practice of approximating a resonant 
cavity with an equivalent lumped circuit in the vicinity of resonance. 

Also, if, in any (or all) terms of (10) or (11), s is replaced by s + 6(6 > 0), the series 
converges as before but now represents a circuit with small losses; K, and w, may, to 
a first approximation, be evaluated as above, and the value of 6 for each term can then 
be obtained from the ‘‘Q” of the circuit at that resonant frequency. 

Another useful expansion is analogous to Foster’s cannonical form for lumped- 
constant reactances. Here it becomes an infinite product. If Z(0) = 0, 








FT = + (s°/2) * 
lint es 1+ eat (12) | 
while if Z(0) = « 
y + 1+ (s°/2 / Qn) 
i= 13 
os Wi + (s°/wz) (13) 


where L or C is the low-frequency inductance or capacitance, Q, are the series-resonant 
(angular) frequencies, and w, the anti-resonant frequencies (both numbered by in- 
creasing magnitude).’® 

The various terms in (12, 13) are much more easily evaluated than those in (11). 
This is somewhat compensated by slower convergence. As a simple example of the 
use of these products in calculation of reactances, consider a shorted line. Suppose we 
had found (say, merely by considering traveling-wave reflection) that the line is series- 
resonant at 2, = nac/l and anti-resonant at w, = (2n + 1)xc/2l. The low-frequency 
inductance is then easily calculated, and substitution in (12) immediately gives the 
reactance at all frequencies. In this simple case, the product can be evaluated for general 
s because it is known to be the infinite product expansion of LZ tanh (ls/c). These same 
principles can obviously be extended to any circuit and will yield impedance formulas 
not easily obtained in any other fashion. 

We turn now to general, non-reactive impedances. The resistance will, in general, 
assume an absolute minimum at some point because Rk > 0. At such a point, the im- 
pedance behaves locally very much like a reactance. The derivative dZ/ds = dX/dw 


is real and positive and even satisfies (9). The lower bound R, of & represents essentially 





183A fter developing (12, 13), the author discovered the recently available article, K. Franz, Das Reac- 
tanztheorem fiir beliebige Hohlréume, Elek. Nach. Tech., Bd 21, Heft 1/2, 8 8, (Jan/Feb 1944). 
Franz proves (12, 13) for cavities of arbitrary shape. For a comparison of his work and that of the 


present author, see Appendix B. 
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a constant series resistance because the function Z — R, again satisfies (8) and the itali- 
cized conditions immediately above. 

Except for an added pure reactance, the resistive part of a generalized impedance 
specifies its reactance completely. If R(w) is the resistive part (on y = 0), the formula is: 


[RQ =RO ay 


“0 i __ Cc 


Xe = = 


This formula was proved for the case of lumped-constant circuits by Bode.’ 

The various other formulas of Bode have not yet been extended for the completely 
general case. The original proofs were based on Cauchy’s Integral Formula and require 
that various integrals over a large semi-circle in the right half-plane approach zero with 
increasing radius. For rational functions, this proof offers no great difficulty, but a general 
impedance will have an essential singularity at infinity. This means that the behavior 
of Z, as s approaches infinity, may vary radically with the direction of approach. 

If, however, Z is bounded on the real frequency axis, then it may be shown that it 
is bounded in the entire right half-plane. This fact immediately makes many™ of Bode’s 
formulas available for this class of functions. 

However, the phase-area law: 


f= 
|, X dilog f) = 5 (R=) — RO) 

is certainly not valid unless R() has meaning. Even if Z is bounded, R() may not 
exist; for example, the impedance of a shorted section of slightly lossy line, tanh (b + s) 
(6 > 0), is bounded on the real frequency axis, but oscillates indefinitely as w — @. 
If Z is bounded and R(@) exists, then Z will approach R() uniformly in the entire 
right half-plane. The phase-area law is then easily proved by Bode’s method”’ for this 
class of functions. 


APPENDIX A 

It is well-known that the characteristic impedances of transmission lines and wave- 
guides may have branch points on the real-frequency axis. We have indicated in Sec. 
2 our reasons for excluding infinite lines or guides, but it might be well to point out 
explicitly why the input impedance of a finite line or guide cannot have such a branch 
point if the terminating impedance has none. 

Recall that the propagation function 6 of a line or guide hasgexactly the same branch 
points as the characteristic impedance Z, . In both cases, the point is such that varying 
the complex frequency once around the branch point merely changes the sign of Z, and 8. 

Now, in the expression for the input impedance, Z, and @ appear only in the combi- 
nations Z, tanh @LZ and (tanh 6L)/Z, . Since the hyperbolic tangent is an odd function 
of its complex variable, these expressions do not change at all when the point in question 
is circled. Therefore, by definition, the input impedance does not have a branch point. 


“Specifically (Bode, op. cit.): 
p 301, No’s I(b), I(d), III(b), IV(b), V(c), VI(b) 
p 335, No’s I, III, IV, V 
Bode, op. cit., p 286. Use (Z — R(@))/s as the integrand instead of Z/s. Formulas II, p 335, are 
also valid for this class of functions. 
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APPENDIX B 


It was mentioned in See. 4 that Franz has proved (12), (13). His method employs 
the theory® of the distribution of the eigen-values of the wave equation. Whén the 
developments of this paper were originally evolved, the author was unaware of these 
theories and so did not use the following condition, which is an immediate corollary 
of the present argument of Sec. 3. 

Lim(n(f)/f) asf approaches ~ exists, where n(f) is the number of resonances below f. 
Franz uses this condition in his proof. Since it may be shown that (8) and the italicized 
conditions immediately above (8) do not imply this condition, the author’s theorem is, 
from a purely mathematical point of view, more general than that of Franz. 

Similarly, the addition of this extra condition does not invalidate any of the present 
results since it merely means working with a subclass of those functions satisfying (8) 
and the italicized conditions above it. It is quite likely, however, that restrictions such 
as that just given could be used to obtain further results not implied by our earlier 
conditions alone. Consideration of slightly lossy cavities immediately suggests that a 
more general restriction may exist on the asymptotic behavior of the number of maxima 


of FR on the real frequency axis. 
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A GENERAL APPROXIMATION METHOD IN THE THEORY OF 
PLATES OF SMALL DEFLECTION* 


BY 
M. ZBIGNIEW KRZYWOBLOCKI (University of Illinois) 


Introduction. The methods developed up to now in the theory of plates do not lend 
themselves readily to the solution of the problem in the case of an arbitrary shape of 
the plate. The papers of Nadai, Navier, Reissner, Southwell, and others, deal only 
with plates of special shapes. In the present paper the author develops a new method 
which aims at the solution of the problem quite independently of the shape of the plate. 
The outstanding feature of this method is that it is entirely independent of the shape 
of the plate; its suecess depends, however, on the type of loading. If there is no lateral 
load, i. e., if the forces act in the plane of the plate, the method is successful independ- 
ently of whether the forces are distributed uniformly or not. Thus, the critical stresses 
of a plate of arbitrary shape can be found. If there is a lateral load, the method is not 
always successful. There are cases in which it is not possible to find a particular solution 
or a fundamental integral of a non-homogeneous linear partial differential equation of 
the fourth order with the aid of the methods now available in mathematics. 

In general, the method involves the following two steps: 

(a) to obtain a function which satisfies the linear partial differential equation of the 
fourth order for the deflection of the plate, but does not satisfy the boundary 
conditions, 

(b) to obtain a solution of this partial differential equation which satisfies the bound- 
ary conditions precisely or approximately. 

In the case of a linear homogeneous partial differential equation of the fourth order 
with constant coefficients, there exists the possibility of separating the variables, and 
a complete set of functions can generally be found with the aid of conventional methods. 
The solution then is a simple series of these functions with coefficients which are de- 
termined from the boundary conditions. In the case of a linear homogeneous partial 
differential equation of the fourth order with variable coefficients, the problem becomes 
much more complicated since, in this case, the variables cannot be separated. A method 
given by S. Bergman [1]f then leads to a complete set of functions, and in particular 
to convergent expansions in terms of these functions. The complete solution of the 
differential equation is given by a double series. 

To keep the possibility of applying the present method to the case where the forces 
acting in the plane of the plate are not distributed uniformly, the author chose Berg- 
man’s method. In order not to obscure the whole procedure by long calculations, the 
forces have been assumed to be uniformly distyjbuted in the plane of the plate in the 
numerical example elaborated in the present paper; i.e., the coefficients in the linear 
homogeneous partial differential equation of the fourth order are assumed constant. 


*Received Aug. 4, 1947. These investigations were carried out when the author was a post-doctorate 
Fellow in the Program of Advanced Instruction and Research in Mechanics at Brown University. The 
author wishes to acknowledge his indebtedness to Dr. S. Bergman, Harvard University, for helpful advice 
and suggestions and to Professor G. H. Handelman, Brown University, for some corrections during the 
preparation of the manuscript. 

+Numbers in the brackets refer to the bibliography at the end of the paper. 
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However, this example illustrates the general method. In the case where the forces in 
the plane of the plate are not distributed uniformly, the general method remains un- 
changed. On the other hand, some changes have to be introduced into the numerical 
procedure and this has been explained throughly in Ref. 10. Thus, with a certain amount 
of additional_numerical work, the present method may be extended to the general case 
of variable coefficients. The boundary conditions are only satisfied approximately with 
the present method. 

A brief resume of Bergman’s method of solving homogeneous linear partial differ- 
ential equations of the fourth order is given in Sec. 1. Since this method represents the 
solution in the form of infinite series, which involve in practical problems many thousands 
of operations of a simple nature, this technique is particularly suitable for punch-card 
machines. Such machines were employed in obtaining the numerical results given in 
this paper. The functions appearing in the infinite series are defined and evaluated, 
and the various partial differential equations of loaded plates with small deflections are 
discussed in Sections 2 and 3. As an illustration of the procedure, the method has been 
applied in Sec. 4 to the calculation of critical stresses in a triangular plate loaded uni- 
formly in its plane parallel to the x-axis. The example shows that once the scheme for 
the calculations has been set up, the computations are quite straightforward and simple 
and may be performed by an inexperienced staff. All details of the example are explained 
together with the method of avoiding solution of determinants of higher order. Several 
final remarks close the paper. 

1. Solution of linear partial differential equations of the fourth order. The proof 
and complete explanation of this particular method of solving linear partial differential 
equations of the fourth order can be found in 8. Bergman’s paper [1]. Only a brief resume 
of this method will be given here. The differential equation considered is the following: 


U,eee2e + MU,, + LU,» + NU, + AU, + DU, + CU = 0. (1) 


The coefficients, M, L, --- C are assumed to be entire functions of z and z* where, 
in terms of complex notation, x + iy = z,x — ty = z*. The function U is also assumed 
to be a function of z and z*. When z and z* are conjugate complex numbers, x and y 
are restricted to real values. If x and y aré allowed to take on complex values, then 
z and z* will be independent complex variables. The asterisk denotes the conjugate 
complex number and partial derivatives are denoted by subscripts. If one sets U = 
u + iv, Eq. (1) is equivalent to a system of two real fourth order partial differential 
equations for two real functions u and v. If for real values of 2 and y the functions L 
and C are real and M = N* = M, + iM,, and A = D* = A, + 7A, , then the first 
of these two equations contains u only and has the form: 


AAu + au,, + 2bu,, + cu,, + du, + eu, + fu = 0, (2) 
— 
where 
ou Ou Ou 
Adu = — + 2-9-54+ 
"Oa! + op oy Oy” 
a= 4L+ 8M,, b = 8M, , c= 4L — 8M,, (3) 


Il 

— 
> 
2 


d = 16A,, e=16A,, f 


1948 PLATES OF SMALL DEFLECTION 33 


In applications, we are mostly interested in real equations of the form (2). However, 
it is simpler to work with the general equation (1) without making any special assump- 
tions about the coefficients. The derivatives employed are the following: 




















U=ut+un, 
.  @U 1fau aU aU 1(aU , .aU 
v= Wai _ WU) py, U1, ev) 
Z 2 \dx ’ ay ‘ U, dz* 2\az * * ay : 
i 2 [1 (20 - 2) a 1/2 (au 4 +2 (UC - 28)I 
oe  @e* L2 \dxr ‘ay /1~ 4Lazx \az oy dy \dx oy 
1/#U ey) 
=i (20 +2), ) 
2rr 2rr 2rr 
v= (zv 9 9 _ £0) 
4 \dx Ox dy = Oy 
1 (aU . PU fu) 
veda (x +25 dy day’ /’ 
Vy _ a (2¥ _#U Uv) 
mee 16 Naz" dx* dy” © ayt / 


S. Bergman [1] has proved the following theorem concerning the solution of Eq. (1) 


[see Theorem 2.1. and Remark 9.1. of his paper]. 

Theorem. There exist four functions E“’"(z, 2*, t), k = I, II, « = 1, 2, which are 
defined for sufficiently small values, say | z| <p, , | 2* | < 2, and for | ¢| < 1, possessing 
the following property: If f,(¢) and g,(¢), « = 1, 2,-are any analytic functions of ¢ defined 


in the neighborhood of the origin, then 


U(z, 2*) = z id {EG #, ord 2(1 — | 
(5) 


+ E'™¢, 2*, dod 4 z*(1 — | — ty” dt, 


is a solution of (1). Conversely, if U(z, z*) is a solution of (1) defined in the neighborhood 
of the origin z = z* = 0, then U can be represented in the form (5) by means of suitably 
chosen functions f, and g, , x = 1, 2. 

This theorem yields a local result; i.e., this theorem states that representation (5) 
for U(z, 2*) is valid only in a sufficiently small neighborhood of the origin. By proving 
a theorem concerning analytical continuation in the complex domain of solutions of the 
differential equation (1) with analytic coefficients, Bergman has shown that representa- 
tion (5) is valid in the large i.e. in a finite domain. He also proved the following theorem 
[Theorem 11.1]. 

Theorem. Let U be a solution of (1) represented in the form (5). Let U be regular 
for (x, y)eR’, BE’ (z, 2*, t) for |t| < 1, (2, yeR* + 4$R**. Then the functions f, , f. , 
9; , Jo are regular in $R*. By }R? we denote the domain obtained from domain R® by 
the transformation z“ = $2. 
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Because of the condition |¢| < 1 in (5) this domain cannot be greater than the 
unit circle. In order to obtain rapid convergence of the series, the radius of the circle 
should not be greater than 0.5 to 0.6 for practical purposes. The functions E“*'” (z, z*, t) 


introduced above have the property that 


E“¢, 0, t) oa E?©, 2, d) ve 1, 


Ez, 0, ) = E°” (0, 2*, )) = 0, 


(6) 
Ex, 0, ) = EY, 2*, ) = 0, 
EQ” (2, 0, 2 = Eo (0, 2*, ) = 1. 
Each of the functions E(z, 2*, t) can be represented in the form 
Ez, z*, t) = P¢, 24) + > '2’P¢z, 2*), (7) 
y=l 


where the functions P”” are to be determined. They are defined by the following system 
of differential equations, which arise from the application of the linear differential operator 
IA(E) and from the solution of the equation L4(Z) = 0: 


D,[P)} — 0, 
DP) = —4D,[P] — 2D,[P), 


1 





(n+2) aie (n) (n+1) 
D,[P ] — n? + Qn 4 3/4 [D, [P27] + (2n + 1)D,[P; ] (8) 
+ DP?) + (n+ })dP™) 
+ NPR. + DPS? + CP™), (n = 0,1, 2---), 
where 
D,(A) = H,.,. + MH and D(A) = LH,. + AH. (9) 


System (8) is a system of differential recurrence formulas for the functions P™. 
After P®, P™, --- P®*” have been determined, P*” can be obtained by solving 
the ordinary differential equation, where only P” and P“*” appear. There exist two 
sequences of functions P“"” (z,; z*) and P“”(z, z*), n = 0, 1, 2 --- , satisfying the 
differential equations (8), n = 0, 1, 2 --- , and such that 

Pz, 0) = 1, PS’ (z, 0) = 0, 
(10) 


P™@,0)=0, P”@,0)=0, m=1,2---), 


and 


0, PY? 0) = 1, 


P ae 0) 


Pz, 0) Pe "(2,0)=0, m=1,2-->). 


I 
—) 
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The conditions and the proof of uniform and absolute convergence of these series have 
been given by Bergman [1]. Two other series for P‘” are obtained by interchanging 
z and z*, by interehanging differentiation with respect to z and z2* everywhere and by 


replacing M and A by N and D respectively, and conversely. 
Let us take into consideration only one function Z. Substituting Eq. (7) into (5), 


one obtains: 


Utz, 2*) = in re 259 air — #)-"? dt 


+f ep G, yy] OS la - 2)" at (11) 


‘ r t'2’7P (z, 259 aif = "> None dt + F 
-1 
Let us assume that the integral 


— iZ | [tz Oa - fy” dt, 


may be expanded into a simple power series of the form 


n=f oes é Ja - ty" at = De. (11a) 
-1 vy=0 
To this end let us assume that the function f(z, ¢) under the integral sign can be expanded 
into a power series of the form 
= Pe ~ ’ _ 
{259 = yx c, “(1 — (11b) 
2 v=0 > 


This series (11b) can then be inserted in Eq. (11a), giving: 


n= ff {5h a - ey” dt = ie f= la ~ )"” at 








/—1 
(12) 
= > 2’. 
v=0 
If we make the change of variable t* = p, then 
d. 
dt = in” (13) 
and 
™ ‘ : oe _ T@ +1/2)P(1/2) 
— 2\ "7-1/2 i a v-1/2 1/2 " 14 
I (1 — #) dt i (l—p) “p ” dp To + D (14) 
Equation (14) can also be written in terms of a Beta-function; namely, 
ry + 1/2)P(1/2) _ ( 1 1) 14 
Tot) ~ A" +taa) eae 
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since B(p, gq) = [T(p)T(q)]/[T(p + q)]. Thus: 





a — oe rv + 1/2)T(1/2) , _ — , 
I, = 2% 3 To +) z= do. ; (14b) 
From Eq. (14b) we find: 
c=Z re + _ (15) 








Tv + 1/2)T(1/2) 


In the first integral J, of Eq. (11) the simple form of the power series (J; = }>?.5 2”) 
will be used, whereas the form (11b) together with the values of the coefficient, found 
from Eq. (15), will be used in the second, third, etc. terms in the series (11). For the 
second term in the series (11), we must consider an expression of the form: 


1 — is aaa 
I, = in Uzf iE = at — vy)” dt. (16) 


Substitute in Eq. (16) for f(z, ¢) the series (11b) together with the values of the coeffi- 
cients, calculated from Eq. (15): 


a+ ro) a es \y : - 
I, = | t'z ere (1 — 4 ed d 
lad | v=0 





- Iv + 1) atl a 2 _ p2yr-172 
-Lietieraa’ J, fa-ora sana 


- rv + 1) _v+i DY + 1/2)0(3/2) 
T(v + 1/2)T(1/2) “ T(v + 2) 








I 
iM 


— + a oF oF 2 (v + 1) 
The third term can be evaluated in a similar manner. We easily obtain: 


f" ee] = Ola - rv + 1) va Tv + 1/2)T6/2) 
2 

















“ tt 8 eae 2)” Ty + 3) an 


v—1 


rs) a - $4 _,, ete. 


~ 94 r(3) v+1v~+2) 40+ DOF 2) 





Of course, the functions P’”’ will be functions of z2* alone, if expressions (12) to (17) 
are functions of z alone and vice versa. Taking all four sequences into account, we find: 


U(z, 2*) = pa es Poy 


=0 v=0 


+ = tPF oer 


n=0 v=0 (18) 


a }2 : RPO a” 


n=O pe 


a] 


+> DK, Poe 
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where 

: 1 1 3 i 

Ky, = 1, Kio = et Koo = 41-2" ; 
ge © S 
Ku = 9°9 Kau = 4°3.3 ; 
= - ae (19) 
Kir = 9°3) Kn = 4°39 ? 
a. oe , 3% 1 = 
A» "sce “~"“i¢tiets ”* 


Certain constant coefficients have to be introduced in the process of integration; these 
should be defined by the boundary conditions. Separation of the real and imaginary 
parts leads to the following form of the general solution of the partial differential equa- 
tion considered: 


U(z, 2*) = ask > > K.P ee” | 


n=0 v=0 


+ alm) > DN KP (ee | 


n=0 v=0 


+ bake ba K.P ere | 
n=0 v=0 


(20) 
+ b,2Im[---] + ¢yRe[---] + ¢y2fm[---] 


+ ad } p K.Prngee| 
n=0 v=0 


+ dam), > > KPum(ger | 


n=0 v=0 


In case Eq. (2) is taken as the starting point of these considerations, one does not need 
to take into account the conjugate sequences P“” since both P” and P“” will give 
the same results. Thus only the sequence P“” need be calculated. 

2. Functions P‘” and P‘’’. Of course, there are many functions which satisfy 
conditions (8) and consequently represent the solution. From Eq. (9), we find that 


D,[P“) = Ps). + MP™ = 0, 
from which it follows that 
P® = cos (M'” 2*) or P® = sin (M'” 2*), (21) 
under the assumption that M = constant. The general solution is 
PP? = A, cos (M"” z*) + A, sin (M'” 2*). (22) 
From Eq. (10) one may easily find that A, = 1 and A, = 0. Thus 


PP“ = cos (M"” z*), (23) 
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and similarly, 


P“® = A, cos (M’” z*) + A, sin (M’” 2*) = re sin (M’” 2*), (24) 
DP] = P2. + MP = —4D,[P] — 2D,[P) 
25 
= —4[Pin,. + MPS] — 2[LPS + AP]. sai 
For P“ = cos (M”? z*) = P™ one obtains 
PO = P%,.=0, PY = —M'" sin (M"" 2), (26) 
Pe. + MP™ = 2L M™ sin (M’” z*) — 2A cos (M"” 2*). (27) 


The homogeneous equation P{?. + MP“ = 0 has the solutions 
P® =u, = cos (M” z*) and P™ =y%= ai sin (M'” z*), 


Any conventional method might be used to solve the non-homogeneous equation 
(27) by finding the particular integral, for example, the method of variation of para- 
meters. One may easily show, however, all the conventional methods have proved to 
be highly impractical and inconvenient for higher terms of the sequences P“™” and 
P“™” . For this reason, the method given by Ince [2] will be used. This proved to be 
shorter beginning with the term P“™. As an illustration, let us apply this technique 
to the solution of Eq. (27). We note that the Wronskian 








Uy U2 
A(u, , Us) = 
ut us 
(28) 
cos (M’” 2*)(M’””)”* sin (M'” 2*) 
aa = 1. 
—M"” sin (M’” z*) cos (M’” 2*) 
The functions V, and V, are given by 
Yy,=- | _ r(z*) dz*, V,= | = r(z*) dz*, (28a) 


where the expression r(z*) represents the right side of the non-homogeneous equation, 
r(2z*) = 2LM"” sin (M'” z*) — 2A cos (M” 2*). (29) 
After the necessary substitution and calculations are performed one obtains 


L 
M'” 





y, sin” (M'”? 2*) — Le* + — cos (M'” z*) sin (M'? z*) +C,, (30) 


V3 


Il 


L sin? (M'” z2*) — Az* — ro cos (M? z*) sin (M'? z*) + C,. (31) 
d 
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Finally the general solution is given by 
y = Viu, + Vu, = P™ = P™™, (32) 


Condition (10) implies that C, = C, = 0; and, finally, 





pe = at sin (M'” 2*) — Lz* cos (M’” z*) — a z* sin (M’”’z*). — (33) 


One can easily show by direct substitution into Eq. (27), that Eq. (33) gives the re- 
quired result. In a similar way, the function P“*” and all the others can be found. 
The table, given below, shows the results. In order to obtain the function P“” one 
has to replace z* by z, M by N, D by A and vice versa. 


P“® : cos (M’” 2*), 


z* sin (M’” 2*), 





rae sin (M’” 2*) — Lz* cos (M'” 2*) — 


1? 


2 2 
per Ab 2*’ sin (M’” z*) + E — Ar |e cos (M'” z*) 


























3M’? 6M 
+laan- an - Sepeore 
+ E= 22 b* cos (M"? 2*) + aa sin (M’” 24), 
er 4 Taz sin (M’” 2*), 
Pum ea con (a? a) — 2 sin (2%) — sane sin (Mi 2°), 
poem | a ~ aay apa fain (0 2) — Sar? on 
. ite 7 Fy as 
7 FA = A,-3 ~ = +. ao ke * cos (M'” z*) 
* Pr ai he +3 a . Jn eal 
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PU" » cos (N'” 2), 


(7711) L ° rl/2 
P : aa Sin (N'” z 


puny , DL sn (yi 2) 4 [z ua coe (N'”? ) 


, ) * eae 
Lz cos (N'”* z) — : -zsin (N 2), 
N 


~~ 











3 N' 6 6N 
D’ L’ 2M N'‘” 2C' asa 
+/ a she ene ai a en <a fesin (NV 2) 
6NN'? 6N'” 3 Fae 
DL 2A oT : oor Az, 
th | 2 ~s k cos (N’** z) + 3 Ne sin (N 2) 


pitt . _ l sin (N'” 2), 


aaa 
21) ) 1/2 , 1/2 4 i ait 
as Be cos (V'*" z) —- — Din (NV 2) - 2 eet -Z); 
N N N’”? NY? 
12 L* D* oe DL , ae 
catia | —— — _ fesi (N'” 2) — - 2° cos (N’’* 2) 
6N'? 6NN'?] _— 





DL 2A pe 
4 | — A | sin (N'”* z) 
3N N'” 3 N'” 


—’ D’ 2M aC aie 
+ EZ — 3 aC. cos (N z) 


*y J 2C ;, 1/2 
“+ | ! + wi - ..s ae megiemeent | sin (N' 2). 
6NN'? 2N’°N'? 3N'? 3NN™” 


Direct substitution into Eq. (18) of the expressions for the sequences P‘"”, P‘’*”” 


pv™, Pp?” taken with a limited number of terms, and subsequent substitution of 
Eq. (18) into (1) show quite easily that all terms with the power z’ or less than 2’, i.e. 
z’~", 2” ete. cancel. The terms with the power z’*' cancel only partly; i.e. some terms 
do not drop out. The higher the power of z’*”, the greater the number of terms which 
do not drop out. A preliminary test of convergence of the functions P'’’ should be 
made at the beginning for values of M, L, C, D, N and A close to those which may be 
used in the given problem. In the case considered one obtains for M = L = C = D= 
N = 1, A = 3 after the necessary calculations are performed: 








z* 0 0.1 | 0.2 0.3 0.6 
pm 1 0.995000 0.980100 | 0.955300 0.825300 
| 
= _ ~_ —|_—__—____—__|______ PAE! Paenen: in | Te 
pa 0 | —0.029640 —0. 116540 —0. 257000 —0. 946860 
Sees Ee Deke Cena mm ace iae Oe 
Pci) | 0 0.100698 0.239742 | 0.367548 | 0.987782 


As one can see, the convergence is fairly good for small values of z*. 
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3. Partial differential equations of loaded plates with small deflections. The partial 
differential equation solved by Bergman [1] is homogeneous, i.e., all forms contain the 
function U’. On the other hand the partial differential equations of loaded plates with 
small deflections are mostly non-homogeneous, i.e., they contain a term without the 
function ’. In these cases one has to perform the transformation of a non-homogeneous 
equation into a homogeneous one, a task which is not always possible with the use of 
mathematical methods available today. Let us consider several cases [3]: 

(a) laterally loaded isotropic plate: 
at... a4, a ‘ 
w+ Qasr + Vat (34) 
Ox x oy Oy D 
where 
q, intensity of a continuously distributed load, 


w, deflection, 


__ Eh 
12(1 — v*)’ 





D= 


E, Young’s modulus, 
h, thickness of the plate, 
v, Poisson’s ratio. 
(b) laterally loaded anisotropic plate: 








a¢ _ 3" ; Aw 
D.— + 2H -57,+ D,—% =¢, where (35) 
ar Ox dy oy 
Eu Eye En Gh’ 
D, =_ 12” D, _— 12’ D, = 12 ; Dz, = 12’ 


H = D, + 2D,,. 


(ce) isotropic plate bent by moments distributed along the edges: 


“4 “4 4 
aw . 0 Ww ow 
3 * tare 5 SG; (36) 
Ox dx Oy oy 
(d) isotropic plate on an elastic foundation: 
“4 “4 “4 
ow ow ow gq — kw ” 
+ 2— 5 + — = ho, where (37) 
Ox" ox oy dy" D 
k = modulus of the foundation (pound per square inch per inch of deflection). 


(e) isotropic plate bent under the,combined action of lateral loads and forces in the 
middle plane of the plate: 


O*w Ow ow j . ow . ow ‘ o’w 
L g—._—.. 4 = ¢ N.>G3 N, = 2N.. = 38 

ax’ v © an" oy oy* D E + Ox ill oy” + * dx dy)’ (38) 

where 

N,, N,, normal forces per unit length of sections of the plate perpendicular to z- and 


se) 


y-directions, respect ively, 
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| shearing force in the direction of y-axis per unit length of section of the plate 
perpendicular to z-axis. 
Setting g = 0 in Eq. (38), we obtain the equation for the deflection of a plate under 
the action of forces in the middle plane only. 
(f) isotropic plate with small initial curvature: 


Ow O*w, dw 1 - O'(we + w,) 
ot + Bartel g + Ne 





ax" dx” dy’ ay* Ox 
(39) 
, F(Wo + ws), oa 97(wWo + “| 
A. Ts eh. Of Ste 
ili oy” FON On ay 
where W, , initial deflection, 
w, , additional deflection. 

In general we may write 

AAw + aw,, + 2bw,, + cw,, + fw = X, (40) 


where X, which is equal in most cases to g/D, may be a constant, or a function of z 
and y, or sometimes is equal to zero. The general procedure is the following: one finds 
any particular solution of equation (40). Let it be denoted by w, . After substituting it 
into (40) and subtracting this new equation from (40) one obtains the homogeneous 


equation, 
AA(w — w,) + a(wr, — Wizz) + 2b(w., — Wisy) 
+ c(wy, — Wiyy) +f(w — w,) = 0. (41) 


Let us consider some special cases: 
(a) assume f = const. and X = const. Put w, = const. = EL; fE = X, w, = X/f, 


(b) assume f = 0, and X, a, b, ¢ are constant. Put w, = 3 const. z’, w,,, = const., 
a-const. = X, or w, = 3(X/a) x’, 

(c) assume a = b = c = f = 0, X = const.; put w, = (1/24) const. x‘, 

4 r 

Ow . X 

—~ = const. = X, or w, = =x". 

Ox 
(d) assume L(w) = X where X is a certain polynomial in z and y with constant co- 


efficients; assume also a, b, c, constant. Thus X.= >> a,,2’y" or X = > a,, cos 
(vz) cos (uy). Let w, be a polynomial in x and y of the same degree as X. Thus 
wo, = A,,2'y" or w, = = A,, cos (vx) cos (uy). Substitution of w, and subsequent 
comparison of the coefficients on both sides yield the values of A,, in terms of a,, . 
In case no such combinations like those given above can be applied, or when the 
coefficients a, b, c, f, X, are not all constant, one must refer to the general theory of 
partial differential equations in order to find any particular solution or any fundamental 
integral. The application of Green’s method may be successful in some cases. The 
particular case AAw = f as well as the method of finding the fundamental integral in 
the case of an arbitrary linear partial differential equation with constant coefficients 
are treated in Reference 4. When the equation is homogeneous, the method of linear 
differential operators explained above can be applied and the deflection of the plate 


presented in the form (20). 
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4. Illustration of the method. Calculation of the Critical Stresses. (a) Formulation 
of the problem and explanation of the procedure. As an illustration, this method will be 
applied to the problem of calculating the critical stresses in an isotropic triangular plate 
subject to a pure uniform compression parallel to the horizontal axis. The triangle is 
assumed to have a right angle and both legs are equal to l. One of the legs coincides 
with the x,-axis and the other with the y,-axis. In this case, Eq. (38) will reduce to 
the form 





d*w o*w O*w ow N, 
ont + ax? oy? + ay! +r ant 0, where int (42) 


Let us transfer the figure into an z, y-system of coordinates such that the length of 
a leg equals 0.6. Keeping in mind that 


62, dw _ dwar _ 0.6 dw = Hw _ (2:8) vw 
. dz, Ox dx, LI dx’ da, ~\I/ dx” _— 
one obtains in the z, y-system the following equation: 


4 q 4 2 2 
dw ° Ow dw l Ow (43) 


art sé” Oy ay* + 0.36 ° dx” 

In case the two sides are not equal, one has to express the length of one side in terms 
of the length of the other, and Eq. (43) will remain the same. To each value of N, = 
const., there corresponds a set of solutions w satisfying Eq. (43); but the boundary 
conditions need not be satisfied. However, the boundary conditions can be written in 
terms of a minimal problem, and the solution w so determined as to satisfy this minimal 
problem. This method is explained in detail below. 

In order to find the first critical stress one has to find the least value from this set 
of w. To this end, we must vary N, and look for that value of N, for which the deflection 
w reaches its minimum value. Suppose now that 7 is multiplied by a number m. In order 
not to change the character of Eq. (43), we shall introduce a new system of coordinates 
v2 OW OW ip Ow Fw. Aw _ dw 


f=a2m”,n=ym”, i. ar m'; ro Doe oe m, a ae m’, ete. 


Eq. (43) will be transformed into 
O*w d*w O*w Paw 
ae 2 oe ee | ee -: Gan a 0. 
ar + 25 ap t ant | 0.36" af (44) 
Thus multiplication of + by m is equivalent to the multiplication of the coordinates 
xz and y by m’” or to the introduction of new coordinates ¢ and n. Assuming clamped 
edges all along the circumference leads to two conditions: 
(i) deflection along the circumference equal to zero, 


$ wdc = 0 or $ w’ do = 0, da = (dt” + dn’), (45) 


ce ce 


(ii) the first derivative of the deflection with respect to the normal along the cir- 


cumference is equal to zero, 
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dw, _ f (sey Ls 
go do = 0, or 7 do = 0. (46) 


e 


If the solution is not very far from the state where (i) and (ii) are equal to zero, both 
conditions can be combined into 


‘ lw \° ; 
L = $ w do + A g (aw) do = min. (47) 


c 


Concerning the value of A, somé authors (References 5, 6, 7) assume that it is con- 
nected with the question of what ratio each of the two conditions w = 0 and dw/dn = 
0 should be taken into account. If the function w satisfies the boundary precisely, both 
w and dw/dn would vanish on the boundary. The particular solution chosen, however, 
does not satisfy the boundary conditions completely. That part of the work which 


arises due to the deviation from w = 0 may be expressed as 
k, = C, ) w’ da (48) 
and that part of the work which arises due to the deviation from dw/dn = 0 may be 


expressed as 
, (dw\’ 
> = ’ — 4 
ke = C, $ (a) da (49) 


where C, and C, depend on the properties of boundary supports. The next condition 


imposed is the requirement that 
oar . Lt (dw\’ . 
k,=k, 16; Cs ) wdse= C9 —] do. (50) 
dn 
Therefore, the coefficient A is 
Am Co/C, « (51) 


This requirement implies that the work of the boundary forces which arises from 
the chosen particular solution divides uniformly into both parts (References 5, 6). For 


the boundary conditions w = 0, d°w/dn” = 0, Bergman cites that on the basis of ap- 
proximate calculations the requirement (50) is fulfilled best if one chooses A = 0.002, 
0.5, or 200. In the case considered in this paper, the value A = 1 was chosen after some 


approximate calculations. 
Equation (47) can also be written in terms of the (x, y)-coordinates in the following 
manner, 


do = (dt? + dn’)'” = (m dx’ + mdy’)'”? = m'” ds, (52) 


L = m'” ) weds + A m'” ) (4) ds = min. 








1948 PLATES OF SMALL DEFLECTION 45 


One more condition must be superimposed on the whole system, namely, that the 
energy of the deflected plate should be positive. In the case considered, it is sufficient 
to assume that the energy 


: 1 ff ,,{dw\ 
E, = > | N fs 2) d¢dn>0O Ol = const., (53) 


Since N, = const. this condition may be expressed as 


x 


\2 ; 2 
E = [/ (<w) dé dn =m If ( o ‘ dx dy > 0 or = const. (54) 
JJ \d& ’ m'”? dx 


It is obvious that in Eqs. (44) to (50), (53) and (54), w is expressed in ¢, 7 system 
of coordinates. In the case considered m = 1, 4, 9, 16, 25 and 36 was chosen. 

The deflection w will be expressed by U(z, z*) in Eq. (20). Since Eq. (42) was chosen 
as the starting point, only the sequences P’” and P“” will be taken into account. 
After all necessary calculations are completed, the imaginary part of U(z, z*) will be 
assumed to be real, since one may easily see from Bergman’s method [1] that for real 
partial differential equations the function uin U = u + iv will be real. The circumference 
of the triangle and its surface will be divided into several parts. In the case of a curvi- 
linearly shaped plate, the circumference has to be transformed to a polygon and the 
precision of the calculation will depend on the number of sides. The coordinates x and 
y are measured for the center of each part or side and the sequences P“” and P“” are 
calculated for a finite number of terms by substitution of these measured values of z 
and y. The unknown coefficients a,, , @,2 , b,, , b.2 , are determined from the minimal 
problem with constraints discussed previously. This leads to a set of n linear equations 





with n unknowns of the form 


aL ak 


- ~—_ = 0 
0a, 0a, 


, (55) 
where \ is the Lagrangian multiplier. These equations will possess solutions a,, and 
b,, different from zero if and only if the determinant of the coefficients vanishes; that is, 


provided 
én + Agi, *°°°°°"° Cin + AGin 
: : : = F(A) = 0. (56) 
Cnt + AGnt ere, ee er Cnn + AGnn 
In order to find the first value of w one does not need to calculate the deflection 
from Eq. (20). It is quite enough to take F(A) into account (Reference 8)*, for consider 


the pair of quadratic forms 


*Chapter 13, p. 169, Theorem 3. 
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Then the function 


@— rv = D> (a;; — AV; 2.2; , (58) 
1 
has as its discriminant, the determinant 
| ay, rr Nb, eee Ain — Abin | 
|: ; : = F(d). (59) 
Ani acid AOni Sees oe, Qnn — ADan 
If X; , Xo , As *** A, denote the roots of the \-equation of the pair of forms ® and y, 


F(X) = 0, then it can be proved that & and y may be reduced by means of a non-singular 
linear transformation to the normal forms 


2 2 . 2 
Ni®1 + Aowe + eeeee Tt Awa » 
(60) 
oe eee + 2. 
Thus we obtain a pair of quadratic forms, one with coefficients \; , Az , As «++ A, , and 


the other with all coefficients equal to 1. Hence we may consider the pair of quadratic 
forms (60) instead of (57). Since one looks for the first value, only 4, has to be taken 
into account for various values of m. The minimum of the function \, = f(m), which 
may be plotted as a diagram, gives the value of m for the critical stress. 

(b) Calculations of the sequences K,,P“™ (z*)z"*” and K,,P°™ (z*)z"*”. Let us take 
the first three terms of the first sequence (18): 








’ ; 1 ort? 
x ie > (710) o* ‘ 4 pil1)) 2 
U(z, 2*), = F (2*)z + 5/ ( eFD 
(61) 
3 >(712) 7x) ie 
uP Sr Ese) 
If we substitute (23), (33) etc. into (61), comparison of (2) and (43) gives: 
ot ékktndasans~et, wd 4, ao we 
“= 0.36 D’ EEE SR, Wee ae eS 
2 T 2 AT 
Lao Ne ye fe ; 
8 0.36 8D 16 0.36 16D (62) 


: rv N,. _ UN, 

16 0.8616D 5.76D’ 

A=D=C=0. 

Let us assume in (1) that M = 1, L = 2, N = 1 and substitute these values into (23), 
(33) ete. and into (61). After necessary calculations have been made we obtain 


‘ ; 2 
U(z, 2*), = (cos 2*)2” sin 2* — z2* cos 2*) ——— 
(63) 
v+2 


I se es ae 
og ee Ct wee 
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47 
Let us assume »v = 0, 1, 2, 3, and decompose (63) into the real and the imaginary parts, 
calling the real part y,, , and the imaginary y,, . Then 





Wo. = xsin x cosh y + y cos xsinh y 


anf 4 2.2 (64) 
+ Ee + <7 —-x—-yt+ Joos x cosh y, 


Yor _ 


—izx cos x sinh y + tysin x cosh y 


‘ ‘ 22 (65) 
+ Be + “— —2-—-y t+ 1 |i sin x sinh y. 

It can be easily seen that y,. can be obtained from y,, changing (sin x) into (—7 cos 2), 
(cos x) into (i sin x), (sinh y) into cosh y) and (cosh y) into (sinh y). Consequently, it 
is enough to give the results for y,, and ®,, only. 

= 0 osl 
6 12 5) 9 + 2 loos cosh y 


4 2.3 5 3 2 
_ ez zy — a Se ss . 
E Sa 2 + y fin 2 sinh y 





























(66) 
2 2 
+ E — W hin z cosh y + zy cos x sinh y, 
6 42 a oe — . 
bn = |= tty — zy 42 y + 2° — y' Joos 2 cosh y 
24 3 
5 5 3.3 3 3 
| > [eta 4 ty 4, 2e'y + Qa? 2zy in ae 
12 6 3 
. = 2 2 Re 3 
+ z= sin x cosh y + sey Y cos x sinh Y; 
7 5. 2 3.4 6 5 3,2 ‘ 4 
» a — ay? — bay’ — 3xy x Ory — day _ | 
Ya = | 40 4 +z 3xy° |cos x cosh y 
, Ba°y rm Baty’ 4 ay? = y’ 
40 
(68) 
4 CS a 
_ bey t zt y y + 3a°7y — v* hin xsinh y 
“ — 62y + y' 
42 x +y 


sin x cosh y + (2*y — zy*)cos x sinh y. 
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In the same way the function U(z, z*), for P ™ can be calculated; namely, 


y+il . v+2 
~ 
ra 


| | 2 Rows 
a »(120) 7)” a piT21) 7m ie pti 2* a r 
UG, 2*)2 = PUME*)e + POM) TT) + 4P MEE sD 


yy 


on Zz" +( 
(vy + 1) (69) 


(sin 2*)z” — z*(sin z*) 


v+2 


» 
te 


; l ag 
+ (sin z* — z* cosz2* + =2* sin :*] arene aoe 
2 (v + 1)(v + 2) 
Let us call the real part #,, , and the imaginary ®,, . Then 


4 1 Ore? 1 B 2,7 
d,, = k t avy +y = + 3y + a x cosh y 


4 2 
(70) 
2 34 
+ ay cos xsinh y — ety - cos x cosh y + “yey sin x sinh y, 
5 3,2 asf 3 
?,, = E + acy + ay _ 2 “J Say. +2 hin x cosh y 
12 
2 , 
+ xy + day’ + yr _ y + y Joos 2 sinh y (71) 
12 
oP ee 
ow 7: Y cos x cosh y += svt at -sin x sinh y, 
) « 
ot+aty—ay'—y dy‘ — 32° a» Es 
o,, = | ——— y a 2 +2 — y |sinz cosh y 
5 9,,3,,3 = 3 ‘ 3 
+ k cds ae yr Sty + 2zy Joos x sinh y (72) 
12 3 
5 cS a 9 — 
% — 43 me Say’ cos x cosh y + - 3x'y so aa = sin x sinh y, 
2 Py? — 5a*y* — 3zy’" Say’ — 2° — 102°y’ 1? 
d,, = k Ty = y ary + sony a poe + 2° — 3xy° |sin x cosh y 


| Bee + ba'y + ay —y 
i - 40 

(73) 
2y — 10a" y” 


32°y — y° |eos x sinh 2 
30 T Ory ‘| y 


f er 4 2 pe 2 4 6 5 5 
a —5xy — dvy +y ZY wy : 
[ — _ . . cos x cosh y + ——— sin x sinh y. 
L 20 | y 5 y 
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c) Graphical representation of the particular solutions. A graphical representation 
for the particular solutions can be obtained in the following manner: 


Let 2 re’. z” = r’ (cos vp + 7 8in vp), 
re ’*. z*” = 7” (cos vp — 7 8iN vd). 


The values of y,, and ®,, can be calculated for some number of points, say for r = 0.2, 
0.4, 0.6, 0.8, 1.0 and @ = 15°, 30°, 45° --- 90°. These points are plotted on a polar dia- 
gram r, @ separately for y,, and ®,, . Points of the equal value of y,, and ,, can be 
obtained by interpolation and the lines connecting these points will give the graphical 
representation of the particular solutions y,, and #,, . Some of the results obtained 
in this way are shown in Fig. 1. Similar graphs can be obtained for P’””’. The real part 
of P” will be the same as that of P’”” and the imaginary part of P“’” will have the 
same magnitude but the opposite sign to that of P””. 

(d) Calculation of derivatives. The derivative dw/dn has to be calculated for Eq. (52): 


dw _ don _ Wen OE | Vm BY _ I | es + 24a (74) 


dn dn ox On dy on 7 91/2 Ox Oy 


since @x/dn = cosa = 1/2”, dy/dn = sina = 1/2”, a = 45°. 
For the horizontal and the vertical sides Eq. (74) changes to dy,,/dn = (dy,,/dy),-. and 
dy,,/dn = (d,,/dx),-o respectively. For Eq. (54) dy,,/dx, (dp,,/dx)z-o and (dy,,/dzx),«0 


must be calculated. The results are given below: 


4 4 2.3 
oven (2° + zy’ — x)cos x cosh y — jet . —-27’- v* hin xz cosh y 
- (75) 
— ysin x sinh y, 
o 4 4 2.2 
oy = [xy + y*° — yjcos x cosh y + k aa 4 . —-2z-y+ 2 os x sinh y 
oy ‘ 
(76) 
+ xsin xsinh y, 
oe = 2 = +ay>—a2+2°y+y° — y)cos x cosh y 
an - 
at+y , vy . his : 
——e + > ~*~ y pine cosh y 
(77) 


4 1 2 n ‘ 
+ (! Ty +22 yy + 2)eos x sinh y 


] 


+ (x — y)sin x sinh ¥ |p 


To obtain the derivatives of Wo, the same changes as explained above have to be in- 
troduced. In a similar way, all the necessary derivatives dy,,/0x, dp,,/dy, etc., Apo, / A, 


Ado, /y, ete. were calculated. 
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Fia. 1. The functions of ¥ and ¢ in polar coordinates (r, ¢). 


Real part — — — — Imaginary part 
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(e) Procedure for calculation—It is easily seen that the number of operations involved 
in the evaluation the above equations may reach a hundred thousand or more. For 
this reason punch-card-machines were used. The circumference and the surface of the 
triangle were divided into some number of parts; for each part, the coordinates of the 
center point were defined and substituted into the above equations for y,, , ,, and 
for all necessary derivatives. 

The calculations were performed simultaneously for six values of m, which is equiv- 
alent to the multiplication of all terms in x and y of the first degree by m’”” of the second 
degree by m, of the fourth degree by m’, etc. Of course, the values of the trigonometric 
and hyperbolic functions also have to be changed with change in m. In the next step, 
the squares in (52) and (53) were calculated. This is equivalent to the squaring of 16 
term polynomials, since 4 groups y,; , ¥v2 , ®,: , ®.. with 4 terms in each group were 
taken into account. The integrations in Eq. (52) and (53) were replaced by multipli- 
cation of each y,, by the corresponding partial length or surface and summation on 
the punch-card-machines. Next, the smallest roots of (56) were found in several cases. 
In order to obtain a good picture of the influence of the order of the determinant on 
the value of \, the value of the smallest root, 4, was calculated for determinants of 
various orders. Determinants were calculated for all six values of m. The results of 
these calculations are given in Table I for three values of m in the neighborhood of 
the minimum. All determinants were calculated by the use of punch-card-machines. 


TABLE 1. , = f(m) 
































Value of 
| Coefficients Order of -——— 
No. Used determ. m = 3? m= 4 m= 5 

1 | ao 1 —54.5 —5.601058 —8.660602 

2 | der, Des 2 2.8 —0. 15674675 —0.962 

3 ao, Ao2 2 | —8.77342465 —0. 163899924 —1.203058734 
yeas: mpgs ates 

{ An, Ao, 4 —0.0394278 —0.000211601 —0.603463788 

ba, Doe 








To obtain more precise value of \,, , the values of the two functions Yo, , Yoo were 


calculated for all integers between 17 and 24. Table I shows that it is sufficient to take 
into account only the first function Yo, for the location of \,, . The results of calculations 
for Yo, are given below. 


TABLE 2. , = f(m 


m 17 18 19 20 

d —5.344245 —5.288552 — 5.225255 —5.178953 

m 21 22 23 24 
—5.306672 —5.329772 — 5.639540 —6.576221 


Thus \,, is at m = 20. To obtain still more precise value of \,, one may repeat the 
calculations for several decimal values between 19 and 21. 
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The value of the critical stress is given by Eq. (62): 


mN = we (mN,) = Neer ; N = 1, 
(78) 
N.cr = 20 X 5.76 7 = 115.2 2 
This may be considered as the first approximation* of N,,, . 

This method may also be applied to the case where the coefficients in Eq. (40) are 
analytic functions in x and y. Since the procedure in this case remains the same as in 
the case of a partial differential equation of elliptic type and second order one may refer 
to Reference 10. 

The whole procedure may be performed by a graphical method. Diagrams for the 
particular solutions, for the first and second derivatives, can be plotted as explained 
above, and the corresponding values read off for each plate shape considered. There 
is a possibility of preparing a set of such diagrams for various values of forces N, , N, , 
N,, , taken in certain intervals, similar to charts used in practice for buckling of columns. 
This technique has been used in order to check the results obtained on the punch-card- 
machines; the agreement was bound to be astonishingly good. 

Similarly, in the case where the coefficients M, L, N --- in (1) are analytic functions 
of z and z*, one must refer to the general method of solution of such equations by use 
of punch-card-machines [10]. 
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*For comparison purposes several values are given below: square plate, built-in edges, compression in 
z-direction, N,., = 105 (D/a?); square plate, simply supported edges, compression in z-direction, 
< , | I - , 


Nicr = 4x°(D/a?); equilateral triangular plate, simply supported edges, uniformly compressed, V,, = 
4x?(D/a?), a = height. 


ON A PROBLEM IN PLANE STRAIN 


BY 
H. J. GREENBERG anp ROHN TRUELL* 


Brown University 


1. Introduction. The present paper concerns the application to a particular problem 
in plane elasticity of the method recently developed by W. Prager and J. L. Synge’ 
for obtaining approximate solutions to boundary-value problems in elasticity. 

The problem considered is that of an infinite bar of rectangular cross-section in a 
state of plane strain caused by compression as shown in Fig. 1. We assume the conditions 
of the test to be such that the bar is compressed by a prescribed amount 2a, the points 
on the top and bottom faces being permitted vertical (u.) but no horizontal (u,) dis- 
placement. Thus, after deformation the cross-section of the bar will appear as in Fig. 2. 

This is not a simple compression test. In fact, the exact determination of the stresses 
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in the bar requires the solution of a mixed boundary value problem with the displace- 
ments prescribed along x, = -th and the stresses prescribed (to vanish) along xz, = +b. 

We propose in the present paper to determine the relation between the force F per 
unit length of the bar, the compression 2a and the Young’s modulus E for the material 
assuming the value of Poisson’s ratio vy to be known. (We note that the range of v for 
most materials is small as compared with the possible range of E. In any case, by assuming 
extreme values of vy one may determine the effect of varying v on the relation between 
E, F and a.) Knowing this relation, the value of Z for a given material can be determined 
from test measurements of F and a. This test has the advantage that the usual precau- 
tions taken to eliminate friction between specimen and end-blocks become unnecessary. 

To solve the problem stated above one needs only to determine the strain energy 
of deformation of the bar. This we shall do to a good approximation by finding upper 
and lower bounds for this quantity. These bounds can be obtained by means of the 
Prager-Synge method. To apply this method we must construct artificial states of stress 
*Received Aug. 30, 1947. This paper is based on a report prepared for Watertown Arsenal under a 
contract in Applied Mechanics. 

IW. Prager and J. L. Synge, Quarterly Appl. Math. 5, 241-271 (1947). 
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which approximate to the natural state of stress. These artificial states of stress will 
be discussed in the following section. 

2. Specification of states of stress. According to the Prager-Synge method a state 
of stress, given by a set of six stress components E,; which are functions of position, is 
looked upon as defining a point or vector S in function space. Given two states S and 
S’, the scalar product of these states is defined to be 


s:$’ = [ Bae: dv 


{; are the strains corresponding to the stresses E/,; of the state S’ and are com- 


where €¢; 
puted from the #{; by means of Hooke’s law; the integral is extended over the volume 


of the body. It is easily verified that 
S-S’=S'-S 


so that the scalar product is commutative. 

The sum S + S’ of two states is defined as the state with the stresses E;; + E/; , 
the product cS of a state by a scalar c is defined to be the state with the stresses cE;, . 
It follows that the usual laws of vector algebra remain unchanged. 

We shall reserve the notation S for the natural state of stress E;; in the bar. Thus, 
the stresses #;; are assumed to satisfy the equilibrium and compatibility equations and 
in addition all of the boundary conditions of the problem. We note that 


S=S-S= [ Buses ae (2.1) 


so that the square of the magnitude of S is equal to twice the strain energy of deformation 
of the bar. 

According to the method, we must now select various artificial states of stress which 
satisfy some but not all of the conditions of the problem. We note first that our boundary 
conditions specify the displacements on the top and bottom of the bar and zero stresses 
along the sides. This makes our problem one of a type considered by Prager and Synge 
and designated by them as Displacement Boundary Condition (DBC). In accordance 
with their treatment of problems of the DBC type we shall need the following states: 

(a) A state S* called the completely associated state. The stresses E*; must satisfy 
the compatibility equation. The displacements defined by the E¥, are required to satisfy 
the prescribed boundary conditions on displacements. 

(b) A sequence of states S{ , Si, --- , Sf called the homogeneous associated states. 
The stresses ,#/; of the state S/ must satisfy the compatibility equations. The displace- 
ments defined by the ,E{,; are required to vanish wherever displacements are prescribed 
by the boundary conditions. 

(c) A sequence of states S{’, S3’, --- , S4’ called the complementary states. The stresses 
2%; of the state S/’ must satisfy the equilibrium equations and in addition the boundary 
conditions on stresses. 

Since the bar under compression is assumed to be in plane strain, the displacements 
u, and uz are functions of x, and z, alone and there is no displacement perpendicular to 
the x,x,-plane. The boundary conditions to be satisfied by the state S* are therefore 


simply 
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u%(z,, +h) = 0 
(2.2) 
us(a, ’ +h) 


where u¥ and u¥ are the displacements which give rise to the stresses E¥, of the state 
S*. Accordingly, the boundary conditions to be satisfied by the displacements ,u{ and 
,uz associated with the stresses ,E{; of the state S/ are 


sui(x, , +h) = 0 


Il 


—a, ud (x, ’ —h) = 4a, 


(2.3) 
puz(a, ’ +h) = 0. 


It is simplest to obtain the associated states by beginning with displacement functions 
satisfying the boundary conditions, computing from these the strains and finally from 
the strains computing the stresses by means of Hooke’s law. The stresses so obtained 
are of course automatically compatible. 
In plane strain the stress components E,, and £,, are identically zero and 
E33 = VW(E, + E 2) (2.4) 


where v is Poisson’s ratio. Since all quantities are independent of the variable x, , the 
equilibrium equations take the form 


;' OB , IE _ 
Ox, ad Of, 0 
(2.5) 


Ox, 0X2 


Since there are no loads applied to the sides of the bar we have the following boundary 
conditions on the stresses 


E,,(4~b, 22) = E,2(+b, x2) = 0. (2.6) 
The stresses ,£%/ of the state S/’ are chosen to satisfy (2.5) and (2.6). Only the three 


ip 
functions ,£{{ , ,.£{i and ,#2; need be chosen for each gq, since ,£3{i = ,E{; and we may 
take .E{f = ,E4 = Oand £33 = v(,Ei{ + ,£23) by virtue of the preceding remarks. 
3. Upper and lower bounds for Young’s modulus. Prager and Synge gave the fol- 
lowing formula for upper and lower bounds on the quantity S’ which is twice the strain 


energy of deformation of the natural state: 
L (S*- 1)" < 8° < Ss” — Do (S*- 1); (3.1) 
q=1 >= 


the states If are obtained from the states S/ by orthonormalizing these latter states 
according to the requirement that 
1, if p=r 
Iy-I, = 
\ 0, if pwr 
and the states I/’ are similarly obtained from the states S{’. Since the strain energy 
of deformation is equal to the work done by the external forces we have 


S? = 2Fa (3.2) 
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or in terms of the average normal stress ¢ = F'/2b exerted on the bar by the end blocks, 


Fs = dabe. (3.3) 


Furthermore, from the definition of the scalar products involved, the conditions on the 
various states S*, S} , S/’ and the manner of orthonormalization of the states S/ and 


S;’ it follows that the upper and lower bounds in (3.1) are proportional to Young’s 
modulus E and to a’. We can therefore write for (3.1), taking into account (3.2), 





L,Ea” < 2Fa < U,,Ea’ (3.4) 
where : 
l ” 2 
L, = 33 > f°"). 
Ea p> ( a) 
(3.5) 
: 1 = ‘ 
ee ee 
et > (S*- 1’) 
From (3.4) we at once obtain 
2F . 2F 
— <£ < — ‘ (3.6) 
aU,, aL,, 
or in terms of the average stress o and the ‘‘average strain” « = a/h 
Ab a <E< 4b o (3.7) 


hU,, € cae AL,, e 


Thus, we obtain upper and lower bounds for Young’s modulus in terms of F and a. 
Alternatively if we are given EF and F we can find upper and lower bounds for a or given 
E and a we can find upper and lower bounds for F. Since L,, increases and U,, decreases 
as more terms are added in (3.5) (i.e., as n and m are increased) the bounds can be 
made to differ by a negligible amount so that in effect we determine the numerical co- 


efficient in the formula 


We note that the quantities L, and U,, depend on Poisson’s ratio » which we must 
therefore assume to be known. The exact form of the dependence however can be easily 
obtained for n = m = 1. 

4. Selecting the functions and computing the bounds.” For convenience we shall 
take b = h = 1 and vy = 3 in the following. The completely associated state S* and 
the five homogeneous associated states S{ , S; , S} , Si , S{ which we use are given in 
Table 1. Tabulated are the displacement components, strains and stresses corresponding 
to each of these states. For convenience we have also set E = 1 anda = 1; these quanti- 
ties are easily reintroduced into the final answers by proper proportionality factors. 

The complementary states S}’, p = 1, 2, --- , are most conveniently defined by means 
of the Airy stress function y. The stresses are then computed from the function y by 


taking e 


2In this work we were ably assisted by Mr. F. Edelman who performed the bulk of the computations. 
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a. ae 

vt Bae | 

° | 

’ Oo { 
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| 

ee = ll ay 

mas ax, dx, } 


When defined in this way, the stresses automatically satisfy the equilibrium equations 
(2.5) for any choice of the function y. We note that not any function y may be used, 
however, since the stresses must satisfy the boundary conditions (2.6) where b = 1. 
The five complementary states we use are given in Table 2. 

In orthonormalizing the sequence of states S{ , --- , S{ to obtain the orthonormal 
sequence If , --- , If the scalar products S/-S/ are required. These are listed in Table 3. 
To illustrate, we have S{-S; + 4.8, Sj-Sj; = S{-S; = 0.925714, ete. Also given in 
Table 3 are the values of the scalar products S*-S/, 7 = 1, --- , 5 and S*-S* which 
are needed in the evaluation of the quantities in (3.5). 


TABLE 1: The associated states. 



























































ins — if | cin let = 1 (Ou out)! ot, = 2 
Jae ‘ ee Psi | —* Qa, 12" 2\az, ' az,/| “ dae 
_Sf |i 2(1 — 23) 0 1 — 23 —2,22 0 
_S | 0 | ~2,(1 — 23) 0 0 33 — 1 
_ Ss fail — 2) 0 3a3(1 — 23) | —ax22 0 

Si i} OO _| —2,(1 — 25)a; 0 —2,2(1 — 2x) xi(322 — 1) 
Ss |iaizsa —2)|_ 0 | Bata — 22) | xixa(1 — 228) 0 

S* jut =0 (I|Iut=-z |eh&=0 i|@&% =O “|= —l 

) es ee 
Y a 3 , | ” 3 , ; a 3 , , 
State Ei, = 4 (2ei, + ei) | : Bis = tis = 422 = q (2629 + ein) 
, | 3 2 | = | 3 
Si 1 | ean wa) a oa ~~ ae s = 4 (1 £2) oie, 
S 3 3, — 1) | 0 3 (37? — 1) 
oa SR = . 
| 
9 > . 9 , 2 
S 9 vill — 2X) | 4 eo | gull — 22) = 
2 oS ee a 
; : oe . | 3 , a cas » 
oe 2(3x2 — |. alt ef = 73(323 — 1 
. 1 eee eee eee a 
ee , | 2s ; | 9 a: , 
Re nil OE eR IO 8 
| 3 | ; 3 
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TaBLE 2: The complementary states. 























































































































oo errs or ay oY 
wae) v Bi= 3 | "ja, ast = oy 
; ¥ 
Si’ | 5 xi 0 0 —1 
| 
? 1 4 2 
2 ~ 42 X 0 0 — TZ; 
ar | 1 2 2\2 2\2 2 2 9..2 
| 9 tal — 24) (1 — 2) 4x,x.(1 — 2%) —2x2(1 — 32}) 
Si’ | 7 (1 —2i) | 8a3(1 — 23)” | 4a,22(1 — 23) —x3(1 — 323) 
Si’ rs 9 tita(l — ai)? | ai(1 — zi)? | —2x.(8at — 4at + 2) | 23(15at — 1227 + 1) 
| 4 ” 4 
State) eff = 5 [2B — Eis) ely = 3 Bi ef = 5 (28s — Bit 
oe. - 8 
sr | = ae 
Si’ | 9 0 9 
2” ; Xi 0 -< ay 
7 8 2\2 2 2 16 2 4 2 2 2\2 
4 2 2\2 4 4+ 4 2 2 
— (622(1 — x1) + 22(1 —= [2x3(1 — 3x;) + 322(1 
. x 4,23(1 — 2) 9 
si’ 3 1“2 1 
— 3z3)] — x)'] 
4 2 2\2 2 4 8 5 4 2 4 2 
g [2ei(l — 2) - 22(152; ~3 22(32; 9 [2a2(152, _ 122; + 1) 
A 
—12z + 1)] —4zxi + x) —2xi(1 — 23)'] 
TABLE 3: Values of S{-S/ and S*-S}. 
ae ae eh | = a eS a a oe 
s; || _3.see6e7 | -08 | 36 | 0 0.537143 | 2 
1] 
_s; | -os | 48 _ = -08 | ~18 0.114286 | 0 
| 3.6 —0.8 L 6.045714 | = 0.925714 | —2. 
a oo 0.8 _| — 865714 
Ss; | 0 = — 365714 | 1.112381 | 0.045714 0 
_ Ss ams 045714 | 0 
S;__||_ 0.537143 | 0.114286 | 0.925714 | 0.045714 .364080 —0.4 
|— a ee eee ee 
s+ || -2 ia | oOo | -04 | 6 
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TABLE 4: Values of ri, = a;,/a,, and of aj, . 
ft gw | ge fF 6g hls | te 
—i=1] 1. | 0.206897 | —0.928571 | —0.190946| 0.007212] 0.258621 
| oe L. 0.011905 | —0.343706 | —0.035183 |_ 0.215774 
ane 1 ow | Ss 0.128711 | —0.164771 | 0.371287 
a ae _ _ L. ~0.044662| 1.940316 
5 | oe a eee - ‘7 4.776450 
TABLE 5: Values of S/’-S/’ and S*-S}’. 
a 2a ae 2 ee a 
si || 3.555556 | 1.185185| 0.948148| _0.948148| 0.135450 4. 
Si | 1.185185 | 0.711111 | —0.496649 | — 0.054180 | —0.045150| 1.333333 
__sy_|| 0.948148 | —0.496649| 8.777141 5.616633 | 0.456424| _- 0 
Si’ | 0.948148 | —0.054180| 5.616633 | _5.238667| 0.247434 0 
. | 0.135450 | —0.045150| 0.456424] _0.247434| _ 1.465237 0 
:. 2 1.333333 0 0 0 6 
TABLE 6: Values of r/i = b;,/b,, and of bj, . 
—___[.p-i[ 2 [8 | «4f]|s fs 
i=1 1 — 0.333333 | —1.123812| 0.113388 | —0.098739 | 0.281250 
2 ~ L. 2.571430 | —0.591654| _0.201460| 3.164056 
a 2 “ = 1. — 0.685644 | —0.039757 |_ 0.155412 
4 ~ a“ - 1. 0.015344 | 0.654785 
oe - i - 1. 0.700065 
Writing 
i = > a;, Sj; p=1 , 5 (4.2) 


where the constants a;, are to be determined according to the usual orthogonalization 


scheme, we have for U,, in (3.5) 


1 


Ve - a 
a 


p=1 


(s*. > a;, § 


In expanding, only products a,,a;, (4, 7 = 1 -:: 


2 
, 
) ‘ m 


1, -:: 


, od. 


(4.3) 


, p where p = 1, --- , 5) of the co- 
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efficients enter. Defining r/, = a;,/a,,, (¢ = 1, --- , p for a,fixed p = 1, --- , 5) we 
have 


on git wt a 
AinA jy — V ipl jpBpp 


, 


where the r/, are obtained for each p as solutions of the (p — 1) linear equations 


p-1 


> (S/-Sir,, = —Si-S., j=l,---,p—] 


i=] 


whose coefficients are the products S/-S/ of Table 3. The values of the ri, are given in 
Table 4 as well as the values of the quantities aj, , azz. , -:~ , a5; . From these values 
and the values of the S*-S/ in Table 3 we may compute U,, --- , U; from (4.3). 

Similarly, Table 5 gives the values of the scalar products of the complementary 
states with one another and with S*. Writing 


I’ = >) b,,.S)’ 


7=1 


we get for (3.5) 





L, = os =. (s*. x b,.8!’) n=1,---,5. (4.4) 


q=1 
The products of the coefficients are given by 


Diadig = TUT Dee 144j=1,°°-,@q q=1,:-:-,5 (4.5) 
where the values of the r// and the b;, are to be found in Table 6. 

Evaluating (4.3) and (4.4) we find (reentering the factors E and aso that the quantity 
Ea’ cancels out) the values of L, , --- , L; and U, , --- , Us given in Table 7 below. 


TABLE 7 


BS we 
m=n=1 || 4.500000 | 4.965516 __ 
2 || 4.500000 | 4.928572 _ 

3__|_ 4.676823 _|_4.920992_ 


(5 places accurate, 6’th estimated) 











. 4 cs __ 4.750435 _ __ 4.890932 _ 
5 4.761603 +.856468 (4 places accurate) 


Entering these values in (3.6) we get successively improved inequalities for EZ as follows 
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TABLE 8 

Iteration dS ; 7 
No. see fe ee 
I 40278 ~ < E < .444444 . 
2 | 405797 = <B< 4444ga © 

| 

| 
| F , a 
3 .406422 - < E < .427641 m 
4 408920 = 4s < 421014 © 
F F 


5 | 411822 < E < .420027 — 


At the first iteration the gap between upper and lower bound is 10.3%. In four further 
iterations this is reduced to 1.99%. Averaging the upper and lower bounds yielded at 
the fifth iteration and retaining only 4 places gives the approximate formula 


E = iso © (4.6) 


where we have established that the coefficient is in error by less than 1%. The corre- 
sponding expression in terms of average stress and strain is 


E = .8318 -. (4.7) 


It is interesting to compare these formulas with those for a test in which motion of 
the bar in the z,-direction is not hindered, i.e. with the case of pure compression in plane 
strain. Labeling the uniform stress and strain by ¢ and € respectively, we at once obtain 
from Hooke’s laws, written for plane strain, the relation 


E=(1—»)*. (4.8) 
€ 
In the case v = } which we have taken, this gives 
E = .888889 ©, (4.9) 
€ . 
O! 
F 
E = 444444 — (4.10) 
a 


where F is the load required to produce the compression 2@ in the simple compression 
test. Comparing these formulas with (4.6) and (4.7) we verify the fact that preventing 
slippage in the z,-direction has the effect of stiffening the bar. 








62 ‘H. J. GREENBERG AND ROHN TRUELL 


Given these formulas we can correct the value of E determined from a test in which 
it is assumed that the bar is in simple compression whereas actually due to frictional 
end effects the bar is in a state of compression akin to that of Fig. 2. Using (4.10) one 
arrives at a value EZ for E which is larger than the true value as given by (4.6). Division 
of these equations gives as the relation between E and E 


E = .9359 E, (4.11) 


so that .9359 is the correction factor to be used. 

It is to be remarked that in practice one can not achieve even an approximate state 
of plane strain in compression unless a very long (in the x; direction) specimen is used. 
However, analogous results to those given here can be obtained for a specimen with 
rotational symmetry. The formulas deduced for this case would be of more immediate 
practical value in compression testing of materials. 


ON THE NATURE OF THE BOUNDARY LAYER NEAR THE LEADING 
EDGE OF A FLAT PLATE* 


BY 
F. CARRIER (Brown University) anp C. C. LIN (Massachusetts Institute of Technology) 


1. Introduction. The customary treatment of boundary layer flow’ furnishes what 
is essentially an asymptotic expansion of the stream function (or an equivalent function 
defining the velocities). This solution is valid, of course, only in a region far downstream 
of the leading edge of the obstacle. Since one knows, however, that to all practical 
purposes a potential flow exists sufficiently far from the obstacle, there remains only 
a closed region in the neighborhood of the origin in which the nature of the flow is not 
known. It is the purpose of the present paper to show the nature of the flow of an in- 
compressible fluid in this region when the obstacle is a flat plate. 

2. General development. The equations governing an incompressible, viscous, fluid 


. . . 2 
flow are, in conventional notation’, 


nga te 1 4 (2) Po yaw, (1) 

uy Se + a 4 (2) 2 = yan, (2) 
uy, Oy _ 

ete oe (3) 


These may be reduced to a single equation in one unknown function in any of several 
ways. For convenience, we first introduce the following changes of variable: p = 
D;/pUo » © = LyUo/v, Y = YrUo/v, U = Uy/Up = OP/dy, v = 0,/U = —Op/dx, x + ty = 
z = re’’. Here, wp is the free stream velocity. With the introduction of these variables, 


Eqs. (1), (2), (3), are feadily reduced to the form 
a*y/ae? a = 5 [ay/az a°y/a2* 5 — ay/de dy /0e 3. (4) 


Here Z is the complex conjugate of z and, we might note now, u + tv = —2idy/dz. 

If we restrict our analysis to the flow past a flat plate such that ast —> —o,y—y, 
we must impose the boundary conditions dy/d2 = 0 when 6 = 0 and when @ = 2z. 
Furthermore, we must anticipate that the line @ = 0 will be a branch line of our solution 
since the higher derivatives of y will be discontinuous along this line. We may construct 
such a function which.is a solution of Eq. (4) in the following manner. We note that 


vo = Age’? +2277 —- 2? — 2) 
5 
Ar cos £ — cos 3¢) ( ) 
"2 2 


*Received Nov. 17, 1947. 
1See, for instance, L. Prandtl, The mechanics of viscous fluids, in vol. 3 of W. F. Durand’s Aerodynamic 


theory, p. 34, Springer, Berlin, 1934. 
2The subscript unity is adopted so that we may use the variable without subscripts for the quantities 


to be introduced later. 
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is a biharmonic function for which @ = 0 is a branch line, along which the velocity 
dy /d2 vanishes. If we write Eq. (4) in the form 
L(y) = L*(y), (6) 
where L is the biharmonic operator and L*(y) is the rest of Eq. (4), then we may define 
L(Y) = 0, 
L(y) sg L* (Yo), 
(7) 
Ln) _ L* (Yo + al s + Wn-1) bean L* (Yo + ste + Wn-2), 


Over any region in which the series 
Y=vtwt-::: (8) 
converges then, we may take y (as defined above) as a solution to Eq. (4).° In this case, 
a particular integral for ¥, is 
:? = A’r*[(sin 20)/40 — (6 cos 0/64], (9) 


and a biharmonic function which when added to ¥;" adjusts the conditions at the plate, 


is 
a? 
y = jog" Lin r(3 sin @ — sin 36) + 0(3 cos 0 — cos 38) 


(10) 


= cnt + am, 


where a + 38 = 1. As we shall see later, yj” « y in the matching region so that the 
choice a = 1 serves our purpose. The function ¥, becomes then 





vi = Ww? + v? = [A’r*/128][6(cos 8 — cos 38) + Inr(3 sin 6 — sin 34)] 
(11) 
+ [A’r’/40][sin 26 — 2sin 6]. 


12 


It is fairly evident that the higher approximations will contain terms in r’’’, r’’, 
r’? In r, --+ , and will also involve coefficients with higher powers of A. As we (again) 
shall see later, A = .083. Thus, it is evident that within some circle this development 
is convergent and represents some flow. It is not immediately evident that the expansion 
for ¥ converges for y < 1, x > 1. However, since our solution is of the Stokes type, it 
is to be anticipated that the leading term of the expansion describes the flow in a suffi- 
ciently thin neighborhood of the plate. Since we do not need to establish such a region 
of convergence, however, we shall not pursue the question further.* 

3. The matching procedure. The solution obtained by Blasius for the velocity profile 
in the downstream region is usually presented in the form 

®This implies that we base the expansion on the Stokes type slow-flow solutions. 

‘We note, however, that the choice of a for Eq. (10) and the corresponding “‘arbitrary”’ coefficients 


in yn , where n > 1, would be governed largely by this convergence consideration. 
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;, 1 ; 
u=f'(n, ve 2a)" [nf’(n) — f(n)), (12) 


where, in our notation, 7 = y/x'’*. In Eqs. (12), f(m) can be given by a power series of 
the form’ 
f = an /2! - a’n® /2-5!+ ---, (13) 


» 


where a = .332. 
Thus, u = an — a’n‘/2-4! = ar” sin 0/(cos 6)'” 


—a’r’ sin* 6/48 cos’ 6+ ---. (14) 


That is 
u = ar’’6 — a’r’6*/48 + --- (15) 


in the neighborhood of 6 = 0. 
The solution of the preceding section implies that the velocity component uw (i.e. 
Re (—2idy/dz) be given by 


u = Ar'”(5 sin 6/2 + sin 36/2) + --- 


Or 


u =~ 4Ar'”6 + yr? + br? In rot + ---. (16) 


Two facts are evident. If these functions are to represent the same phenomenon in 
any overlapping regions where r is not large, then 4A = a. Furthermore, we note that 
these two functions can never be identically the same since the r’ In r term is lost in 
the asymptotic representation. The latter fact is quite in order since neither solution 
has the pretense of being a complete solution. Therefore, we may establish validity of 
the leading term of this representation (i.e. ¥)) if we can find a region in which both 
and f are valid approximations to the rigorous solution. That is, we must find a region 
in which L(y) — L*(W) « d*y/dr* (say). Our method of solution and its results indicate 
that r ~ 1, 6 < 1 is such a region. But for @ sufficiently small and r of order unity, 
x'’*f(n) is almost precisely the same function as y, (see Eq. (15)) and hence is as valid 
a solution. This, of course, can be worked out by computing L[2z'f()] — L*[2'f(n)], 
but this is really not necessary. Hence, we have that in the neighborhood of r = 0, 
Eq. (5) demonstrates the nature of the flow field provided A = .083. Figure 1 shows 
the four regions into which the flow field may be divided and the solution applicable 
to each. In region IV, it is doubtful if any reasonably simple analytic expression could 
be found for y. However, it is suggested that the flow field could be computed if one 
combined some intelligent ‘‘first guessing”’ (interpolation essentially) with the relaxation 
procedure and applied these to the region IV. 

4. A modified Blasius solution. The alternative procedure’ for developing the fore- 
going solution leads to a slightly modified boundary layer solution which establishes 
more clearly the overlap of the slow-fl6w solution (which includes the leading edge) 
and the boundary layer solution. 


‘Actually, the two authors arrived at the solution independently by the two procedures described 
here, but felt that a joint presentation provided the most lucid account of the matter. 
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(I) Potential Solution (II) Blasius Solution (III) Slow Flow Solution 


It seems that the most natural coordinate system to use for the flow past a flat plate 
is the parabolic coordinate system (é, ») defined by 


x + iy = (& + in)’ (17) 
or 
& + in = r'” exp (0/2) 


By using this system of coordinates, it is possible to give a unified discussion of the 
boundary layer solution and the slow motion solution in the neighborhood of the plate, 
including the leading edge. 

In the parabolic system of coordinates, the equation for the stream function y becomes 


Lip) = (& + 1°) (Wed, — vy Ave) — 2(mb_ — Evy) AY, (18) 


where 


Lip) = (& + Ady — 4(Edy, + nay, — Ay). (19) 


From a consideration of the diffusive effect of the viscous forces, it is obvious that the 
effect of viscosity is essentially limited to regions with 7 ~ 1, close to the solid boundary. 
For very small values of n, the motion is slow. Since y = 0 and dy/dn = 0 at n = O it 
is obvious that in this region a good approximation to (18) is 


L(y) = 0. (20) 
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A polynomial solution satisfying the boundary conditions y = y, = 0 at » = Ois 
vy = Aén’ 
which agrees with Eq. (5) of Sec. 2. This solution can also be obtained from Eq. (19) 


by writing ¥ = n° (f(€) + O(m)) for small 7. 
If & is large while » ~ 1, one may introduce a new variable 


v=-< 


where ¢ is a small parameter.° In the limit « — 0, the equation (3) becomes 


(21) 


(22) 


gE = YPM — WUE) + WOE (28 
] 


(24) 


where 


yo “a ey 


This equation may be solved by putting 


y = ef (n); (25) 


then f() satisfies 
(26) 


a = aaa a }" 


or 
f= 59"; (27) 
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(I) Potential Solution (II) Modified Blasius Solution (III) Slow Flow Solution 








°To give e an explicit meaning, we consider the region » < é. 
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this is the usual boundary layer equation except that 7 is not quite the conventional 

parameter.’ Higher approximations may be obtained by solving (18) in powers of e. 
Thus, the regions of approximate validity of these solutions are as follows: 


(7) potential solution: n> 1; 
(iz) boundary layer equation: 7 ~ 1, » < &; 
(iit) Stokes equation: n <1, irrespective of &. 


These regions are indicated in Fig. 2. It is obvious that there is an overlap of regions 
(iz) and (2z) where the solution (5) agrees with the boundary layer solution. Again, we 
have not established the convergence of the slow flow solution. However, the physical 
considerations which lead to the slow flow and the boundary layer solutions indicate 
that they are at least valid asymptotic approximations, if not convergent processes. 
This is further strengthened by the identical nature of the leading terms of the two 
types of solutions for 7 < 1. 


"Here, 7 = 7/2 sin 6/2, conventionally 7 = r!/? sin 6/(cos @)'”2, 
/ ’ 2 





THE APPROXIMATE NUMERICAL SOLUTION OF THE NON-HOMOGENEOUS 
LINEAR FREDHOLM INTEGRAL EQUATION BY RELAXATION METHODS* 


BY 
F. S. SHAW 


Brown University 


Summary. Relaxation methods are applied to the problem of finding an approximate 
numerical solution to the non-homogeneous linear Fredholm integral equation. As an 
illustration of the technique the deflection of a simply supported single span beam 
subjected to both normal and end loads is found. 

1. The approximation to the integral equation. Consider the equation 


e(2) = fla) +2 f K(x, ded at, (1) 


y(x) being the unknown function. 
Let the interval [a, b], of length L, be subdivided equally into n parts as in Fig. 1., 
so that each subinterval is of length L/n = I say. 





t——}—} 4 
Q=a a 2 : n-1 b 





b -o 


4 | £ | 
= 
or. 
Fig. 1. 


Then equation (1) may be replaced by the approximation 


eh) = fih) + TE {K(h, a)e(a) + K(h, b)p(b)} + > Kh, awa, (2) 


where g(h) represents the value of g(x) at the point x = h, h = a, 1, 2, --- , b, and 
the integral {;*' A(x) dx has been replaced by the approximation 31(A(h) + A(h + 1)). 
To obtain an approximate solution to equation (1) it is necessary, then, to solve the 
set of n + 1 simultaneous equations (2). This may conveniently be done by using the 
Relaxation technique.** 

2. The residual and the relaxation operators. As a preliminary to solving, approxi- 
mately, equations (2) it is necessary to construct two operators which, in the language 
of relaxation methods, are known as the residual operator and the relaxation operator 
respectively. 

Let it be supposed that by means of an inspired guess (or by some other means) 
approximate values have been obtained for all the g(h). In general, the values will not 
be correct, i.e. will not satisfy equations (2), so that the error at each point is given by 


*Received Dec. 17, 1947. ; 
**Some knowledge of the technique is assumed, and it will not be discussed here. For further infor- 
mation see References 1, 2. 
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F(h) = f(h) + al 3 {K(h, a)e(a) + K(h, b)g(b)} + > Kh, Aeli)| — g(h), (3) 


F(h) being known as the residual at the point h,h = a, 1, --- ,n — 1, b. 

Obviously, for an approximate solution of equations (2) we require that each F(h) 
be (approximately) zero. To achieve this use is made of a relaxation operator. 

Consider the effect of making an alteration of amount Ag(3) to the value of ¢(3) 
say. The resultant change in the residual at point 3, i.e. F(3), is given by 


OF (3) 
d¢(3) 





AF(3) = Ag(3), 


[ALK(3, 3) — 1A¢(B) 
from equation (3). At the same time, and as a consequence of Ag(3), there is also a 
change in F(2) of amount 


AF(2) = re Ae(3) 


ALK (2, 3)Ag(3), 


with similar changes taking place in the residuals at the other h stations. 
More generally, the following effects result: 


AF (h, j) = [ALK(h, j) — 6x;J4e(Q) 


| 
AF(h, a) = 4A1K(h, a)Ag(a) | 
AF(h, b) = 3A1K(h, b)Ag(b) | 
AF(a, a) = [4A1K(a, a) — 1]Ag(a) | 


AF(a, h) = ALK(a, h)Ag(h) , fh,j = 1,22, --- ,n — 1, (4) 


AF(a, b) = 3\1K(a, b)Ag(b) 


| 
| 
| 
AF (b, b) = [3A1K(b, b) — 1]Ag(b) , | 
AF(b, h) = A1K(b, h)Ag(h) | 


AF(b, a) = 3\1K(b, a)Ag(a) ~ 


where AF(h, 7) is the change in the residual at station h resulting from a change in 
the value of ¢, of amount Ag, at station j. 6,; is the kronecker delta*. In effect expres- 
sions (4) may be regarded as operators, and are known as the relaxation operators. 
For any arbitrary alteration to the value of the unknown function ¢ at any one station 


(lifh =, 
— { 
(0, ifh ¥ j, 


agi 


- 


ree 
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they give the resulting alterations to the errors in the solution at all stations. Thus they 
provide a means whereby the errors, or residuals, may be systematically’ reduced to as 
small a value as desired. 

For actual use these operators are evaluated in tabular form. On putting Ag(h) = 
1, and dropping the A from AF(h, 7) we have the unit operations Table I. 


TABLE 1. RELAXATION OPERATORS 























Ag(h - 1, : = Alteration to Residual F(/). 
all the rest ——— | a 

zero. F(a) | F(1) F(2) | oe Fin — 1) F(b) 
h=a ALK (a, a) — ) = 1[puxa, a) |ZALK (2, a) 5 a ° AIK(n — 1, a)|3NK(b, a) 
—— oe a eer SN 
h=1 LK (a, 1) | NKA, NK (2,1) eee | AK (n — 1, 1) ALK (b, 1) 

ee |- a eee eee) A se ee 

h=2 ALK (a, 2) NK (1, | NK (2 2) - il, tee | NK(n — 1, 2)) ALK(b, 2) 











h =n —1| NK(a,n — 1)| NK(1,n — | NIK (2, n — »} 
n — 1) — 1) NK(b,n — 1) 











(Z\lK(n — 1, ,b) 4xLK(b, 6) — 1 








h =b M LK (a, 5) |ZALK (1, b) ‘|lpIKe, b) 


3. An example. In order to illustrate the method an example will be worked. Consider 
a single span beam simply supported at the ends and subjected to both normal and 
end loads. Both the normal loads and the moment of inertia of the beam may vary in 
any manner whatever. It is required to calculate the deflection of the beam. 


w(x) 











B 


my 


With the notation shown in Fig. 2., and using the Bernoulli-Euler theory of bending, 
the governing differential equation* is 
*Fig. 2 is diagrammatic only. Actually it is assumed that the neutral plane of the beam is straight, 


that the beam cross section is symmetrical about a vertical axis, and that the force P acts collinearly 
with the resulting straight line axis. 
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dy M(x) 
dx I(x) 


which, for the simply supported beam shown, becomes 


d’y l 
x *y + M,(ax)). (5 
dx Ii I(x) (Py + M,(x)) ») 


M(x) is the bending moment associated with the normal loading system w(x), and so 
for any given w(x) may be considered as known. /(x) is the relevant moment of inertia 
of the beam cross section and may also be considered as known. 

The integral equation appropriate to the problem is 


ol 


| K(x, PAQM A) dt + | K(x, HPACYE) de, (6) 


y(x) 


where A(£) 1/E I(é), 
1é(1 a/L) for O « é <. g. 
K(a, &) 
i ar(l §/L) forza <§ < L, 


and M,(a#) can be evaluated for each point x once w(x) is specified. 
The integral equation is, therefore, of the form 
ol 
y(x) f(x) + A | K(x, &)y(é) dé 
Jo 
where f(x), \ may be considered as known. 
For comparison of accuracy it is convenient to choose a particular problem for 


which an exact solution is known. Accordingly it will be assumed that the load w(z) is 


a uniformly distributed load, and that the moment of inertia is constant. For this 

particular case equation (6) becomes 

w Cor 
\ 4 se 9 4 , ‘s ¢ food 
y(t) = ae (a + L's — 212°) 4 Tal K(x, t)y(t) dé, (7) 
K(x, &) being as given earlier. 

As with most numerical procedures it is desirable to make the quantities concerned 
non-dimensional. To do this, put 
2,L, y = y,L, PL*/EI = A, wh*/EI = B. 


On dropping the subscript 1 (for simplicity) equation (7) becomes 


y(a) B re (a +1— 22°)+ A | K(a, &)y(&) dé, <2 <4. (8) 


he ot 


For convenience* we choose \ = A 2.5, B = 0.5, and take n = 10 so that 1 = 1/10. 


The relaxation operator table may then be completed without further discussion, the 
result being Table IT. 


*The exact value of the deflection at the center of the span is given by 


~ 3 2 
y 5 wL* 24(secu — 1 — fw) 
L 384 HI 5u* : 
where u 2,secu— o. Thus for u = w/4 say, PL*/El 


= (L/2)(P/E1)'. Asu — 0, secu — 1, asu — f, 
= A = 77/4 = 2.5. Also, with y/L = (5/384)wL'/E] = 


< 


1/150, B = 1/2. 
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To compute the initial residuals, note that if the initial solution is taken to be y = 
0 everywhere, then F(h) = f(h), h = a, 1, --+ , 9, b. This gives the first line in the re- 
laxation table which follows—Table III. The relaxation process is then carried out in 
the usual manner. It will be seen that, since the operators are extremely powerful, 
convergence is very rapid. 

From the nature of the problem solved here, and by inspection of the relaxation 
table, it is obvious that advantage could be taken of symmetry to cut down the work. 
It would merely have entailed combining operators 1 and 9, 2 and 8, etc., in Table II, 
and then only relaxing the residuals F(a) to F(5). For clarity however it was thought 
advisable not to do this. 

4. Accuracy of solution. The exact value of the deflection is given by’ 


y_ wl’ | sos — 2x/L) _ | wh 4 
L~ 16 EIu 008 w 1) > ser 





x), 


where u = (L/2)(P/EI)™”. 
Making this expression non-dimensional by use of the previous substitutions, and 
evaluating it for points h = a, 1, --- , 5, the following results are obtained: 

















——— = h=a | 1 2 | 3 4 5 (center) 
———$—<—_— —$<— | —_—_—_|__  URERREE GEE 
Exact Solution. 0 | 0.002,732 | 0.005,175 | 0.007,094 | 0.008,308 | 0.008,727 
Approximate Solution. 0 0.002,737 | 0.005,185 | 0.007,106 | 0.008,331 | 0.008,751 


The agreement is seen to be very good. 

If greater accuracy were desired this could readily be obtained without increasing 
the number of subintervals by using a better approximation for the integration process. 
Thus, for example, Simpson’s three point integration formula* could be used at the small 
expense of increasing slightly the initial labour necessary to compute the relaxation 
operators. However for some problems, depending on the terminal conditions, it might 
be desirable to decrease the lengths of the subintervals adjacent to the two ends a and b. 

It is evident that the method as given here may be applied equally well to obtain an 
approximate solution to the non-homogeneous linear Volterra integral equation. 


REFERENCES 
1. R. V. Southwell, Relaxation methods in engineering science, Oxford Univ. Press, 1940. 
2. F. 8. Shaw, An introduction to relaxation methods (approximate methods of numerical computation), 
C. 8. I. R. Div. of Aero. Report 8 M 78, Sept. 1946. 
3. S. Timoshenko, Theory of elastic stability, McGraw-Hill Book Co., Inc., 1936. p. 8. 
4, R. V. Southwell, Relaxation methods in theoretical physics, Oxford Univ. Press, 1946. p. 15. 


*The use of more complicated integration formulae do not always lead to greater accuracy than that 
obtained by use of the simple trapezoidal rule. See, for example, ref. 4. 
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—NOTES— 


ON BERNOULLI’S METHOD FOR SOLVING ALGEBRAIC EQUATIONS* 
By BERNARD DIMSDALE (Ballistic Research Laboratories, Aberdeen Proving Ground, Md.) 


1. Introduction. Both Graeffe’s and Bernoulli’s method for the solution of algebraic 
equations have the considerable advantage that no initial approximation to a root is 
required. Graeffe’s method, although the convergence is quadratic, leads to factors’ 
whose roots are large powers of the roots of the original equation, and for which the 
determination of angle is therefore quite difficult. Bernoulli’s method, as amplified 
below, leads to factors of the original equation and therefore avoids this difficulty. While 
it is true that the convergence is slower (doubling the number of steps is roughly equiva- 
lent to one step of the Graeffe process) this is not a consideration of first importance if 
high speed computing machinery is available. Also, there is the possibility of applying 
rapidly convergent processes when a “sufficiently good” approximation is obtained. 
Furthermore, the method applies equally well to equations with complex coefficients. 

Bernoulli’s method consists of defining a sequence s, in terms of the coefficients of 
the original equation, and taking quotients s,,,/s, . If these quotients converge, they 
converge to the root of largest modulus. If not, there are at least two roots of largest 
modulus. It will be seen that in this case successive systems of linear equations can be 
defined in terms of the s, whose solutions will tend under certain conditions to the co- 
efficients of a factor of the original equation. 

2. Preliminary remarks. Let 


N 
N-» 
Pre) = Daz”, a=! (1) 
v=0 
with unequal roots z, , 22, --: , Zw: , of multiplicities c, , c. , «++ , Cw , respectively. 


The notation is chosen so that no modulus exceeds any preceding modulus. Let 


N 
Py(z) = Diaz’, a=1, M<N’ (2) 
vy=0 
with roots z, ,2., °°: ,24,andM < N’ if zy. = 0. Let 


N’ M 
Sh - yx Cie ’ 8, —s =. 2) 
y=l v=l 
for all’ n. We then have Newton’s identities for s,, : 


& = N, 


n—1 
> a,s,-, + na, = 0, nmil,---,N—Il, 


v=0 


N 


) » as.., =D, n>N. 


r=0 


*Received Aug. 8, 1947. 

\Bodewig, E. On Graeffe’s method for solving algebraic equations, Q. Appl. Math. 4, 177-190 (1946). 
The methods described there apply only to equations with real coefficients. 

If zy, = O, then so is not defined. Setting s) = N, as we are about to do, however, leads to no error 


in the rest of the sequence. 
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by which s, can be calculated directly from the coefficients of (1) for all* n. The same 
identities hold for s,- when N, a, and s are replaced throughout by M, a, s’. 


Consider the system of equations linear in a,, 


M 
ys Qyn8y>-, = O, a, = I, (3) 
where p = n,n + 1, --:,n + M — 1;n > M. It will be shown in the next section 


that if | zy| > | z4i|, then the system has a unique solution for all sufficiently large 
n and that a,, — a, = O(244;/Z)", 80 that a,, — a, asn— o. Hence for all such values 
of M a factor of (1) in the form (2) is obtained, and (1) can be separated into factors 
each of which has roots of equal modulus. Then for each of these factors a shift of origin 
will result in further factorization. 

3. Convergence. Let 








8n-1 8n-2 a Sd 8n-M 
Aon(8) = Sy 8n-1 a 8n-—-M+1 
Sn+M-2 8n+M-3 wevies Sn-1 
and let A,,(s) for v = 1, 2, --- , M be the same determinant with the column headed 
by s,-, replaced by 8, , 8n41, °** » 8n+a-1 , SO that if Ao, (s)# 0, then a,, = —A,,(8)/Ao,(8). 
Let : 
1 1 ey 1 
B,(Zx1 » 22 °° y Znae) = | 2r1 Zy2 vee Zxm |; y=0,1,---, mM, 
M M M 
2x1 2r2 ah 2M 
where the row with z”~’ is absent. For brevity let 5, = 6,(z,; , 2, +: , Za). We note 


that 55 is the discriminant of (2) and is therefore not zero, since all the roots of (2) are 
different. Then, by factoring, 


n—-M 
A,,(8’) = + (2, ’ 22 “Die Zm) 605, (4) 
and, using the theorem for addition of determinants 


A,,(8s) = + Zz, (CyrCre *** Crm)(Zar °° Zam)” So(2n yt? yam) 


6,(Za1 ’ ideals » Zam) 
= (€;C2 ti Cy)A,,(8’) = = } 


where >, denotes summation over all combinations of the M integers \, < A: < +++ < 
\w chosen from 1, 2, --- , N’, and >.’ is the same sum with (1, 2, --- , M) omitted. 
From (4), we have Ap,(s’) # 0 for all allowed values of M. Since the z’s have been 


?The original Bernoulli method permits an arbitrary selection of the first N of the s, . In case there 
are two equal roots of largest modulus and the remaining roots have smaller moduli, it can be shown that 
for an arbitrary selection the convergence of 8,4:/s, will in general be O(1/n) whereas for this choice it 


will be in our present notation, O(z2/z)". 
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arranged in order of non-increasing moduli, it follows that if |zu.i1| < |za|, then 


Aon(8)/Aon(8’) = (e:€2 +++ Cu) + O(2u4;:/2m)", and therefore that for sufficiently large 
n, Aon(s) # 0. Hence, for such n, a unique solution of (3) exists. If 6, is not zero, we 
can also write A,,(s)/A,,(8’) = (c,¢2 «++ Ca) + O(274:/24)" from which 


A,n(8) /Aon(8) _ {1 + O(2441/2m)")Ayn(8’) / on(8’), 
that is 
a, = [1 a O(2+1/2mu)")a, ° 


Therefore, a,, — a, = O(zy4;/2m)". This equation is easily established if 6, , that is, 
a, , 18 zero. 

It might be remarked finally that if M > N’, the system (3) is dependent for every 
n;if M = N’, then a,, = a, for all n if zy. ¥ 0, and if zy. = O then in this case (3) is 
again dependent for all n. 

4. Example. We apply the method to the equation 


2° — 3c* — (2+ ae* + (12 + 5a)2”’ — (8 + &i)z2 + 44 = O, 
whose roots are 
z, = —1— (1+ 7)” = —2.098684113 — 0.4550898608:, 
2 = 2, z= 1, 
—1+ (1+ 7)'” = .0986841134 + .45508986087, 
and ¢, = ¢; = ¢, = 1, c. = 2. Here | z, | = 2.147459380, | z,| = .4656665498. We ob- 


2 
Il 


serve that | z,|/| z,| = .931, | 2; |/| z.| = .5 and so we would expect the quotients to 
converge much more slowly than the solutions of the second order system. We observe 
that as a matter of fact convergence of the quotient is not visible for n < 30, but that 
n = 30 gives us nine place accuracy in the second order system. 


On solving the quadratic 
Z + a; ,202 + a2,29 = 0, 


we obtain 


z, = —2.098684113 — 0.4550898632, Z. = 1.999999998 + 0.000000004:. 


The third order system is not computed, because significant digits are lost in the 
computation. If we carry ten significant figures throughout, we will have only one or 
two left toward the end of our computation for n = 30. This corresponds to the fact 
that the second order system has converged to eight or nine significant figures by this 
time. This is true in general of the behaviour of these systems, and therefore one would, 
in practice, only carry the method to the lowest order which converges with sufficient 
rapidity (“‘sufficient” in this case referring to the speed with which these computations 
can be performed). On removal of the factor obtained, one could repeat the operation 
with the reduced equation. It is interesting to observe in our present example that 
since the quotients converge, we will, if we go far enough, begin to lose significant digits 
in the second order system, but that in the meantime we will have obtained all the 
accuracy possible with ten digits. 

For any n, the fourth order system has the solution a,,, = —1, a,, = —4 — 1, 
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Qs t + 37, a, = —2i, and the corresponding equation has the roots z, , Z2 , Zs 2 . 
T hie occurs because of the presence of the double root z. = 2. For the same reason, all 
fifth order systems will be dependent. 

A practical difficulty may arise with fixed decimal point machinery in that s, may 
go off scale either on the right or the left. This will happen for maprpeed large n when 


lz, | + This difficulty may be overcome by a transformation z’ = z/r, where r = 
lay. |", M being the order of a convergent process, n being “sufficiently” large. In 
our example we have r = | a», |'””. For n as small as three, | s, | would remain in a very 


satisfactory region. 


AERODYNAMIC FORCES ON A SLOTTED FLAT PLATE* 
By EDMUND PINNEY (University of California, Berkeley) 


The derivation of the flat plate approximation to the aerodynamic forces exerted 
by an incompressible fluid on an airfoil is simple and is to be found in most textbooks. 
The present note treats the less simple case in which the airfoil contains slots oriented 
parallel to its length. This problem is of importance in the theory of suspension bridges 
where slotted roadbeds are used to cut down aerodynamically driven oscillations. 

The spaces between the slots will be called “lanes”. Let V be the velocity of flow, 
and a the angle of attack, assumed to be small. Let the plate occupy the interval —b < 
x <b, y = 0 of a rectangular coordinate system, and let u(z), v(z) be the velocity com- 
ponents of the fluid parallel to the x and y axes, sempeaiivdle, at the point z = x + iy. 
Then, if powers of @ higher than the first are neglected, 


u(z) > V, v(z) — Va as|z|— o. (1) 


By two-dimensional potential theory, v(z) + iu(z) is an analytic function of z away 
from the plate. Consider 


N 
o(2) + iu(e) = iV + Va TJ, @ — tb)'*@ — 1b)” (2) 
1 
where z = t,b and z = J,b are the trailing and leading edges, respectively, of the n-th 


lane, and N is the total number of lanes. The lanes are numbered from left to right so 


that 
l= -1, ae se ces Oe ey ioe we. ty = 1. (3) 


As z— o, it is seen that (1) is satisfied. 
In the n-th lane 1,b < x < t,b, and by (2), 


N 
v(x + 10) + tule + 0) = iV + iVa [], | 2 — tb |'? | 2 — 1b |-”. 
1 
Therefore 
is 


u(c + 100) = V+ Va [J], | 2 — tb |'? | 2 — 1b |, v(x +i0)=0. (4) 


*Received Sept. 22, 1947. 
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Therefore the velocity distribution of (2) satisfies the boundary conditions on the lanes. 
Also the velocity is finite when z = ¢,b, so the Kutta condition is satisfied on the trailing 
edge of each lane. 

By Bernoulli’s theorem, 


p+ 5 ou’ + v*) = const., 


where p is the pressure and p is the density of the air. The quantities p and wu are dis- 
continuous on the lanes. Let p, , uw, be their values on the upper side, and p_ , u_ their 
values on the lower side. Then 


Dp. —p-+ : p(u. — u2) = 0. 


By (4), on the lanes, 


N 


Dp. — p- + 2pV'a JI, | 2 — tb |’? | 2 —1,b[-'? = 0. (5) 
1 
On the slots between the lanes, p, — p_ = 0, of course. 
Let 
N 
&(z) = [|], (¢ — t,b)'(e — 1,b)"'7/(z — b), (6) 
1 
N 
¢(z) = — J], | 2 — tb |’? | x — 1b |'7/(b — 2), for x on a lane, 
(7) 
= Q, for x on a slot.* 
Then, anywhere on the.roadbed, 
Pp. — p- = 2pV'*alb — x)d(z). (8) 


We may now calculate the aerodynamic lift L and moment M from 


b 6 
L= / (p. — p-_) dz, M = (p. — p-)x dz. 
—b v—b 


. 
It is possible to obtain these and higher aerodynamic moments, if desired, by evaluating 
I eS ~ , b . > 


J (@. - v/a = re/6)) ae 


and expanding the result in powers of r. 
If0 <r < 1, by (8), 
} ab 3 ” ( 1 } ab b 
= | P+ — P- dy = Qnpl *ba (1 — 1) a | $(x) = = $(x) az} 
oF ily OT 5 


rJ,z—b/r x — b/r 
NowT 


*Pinney, E., Annals of Math. 47, 221-232 (1946). See § 3. 
tPinney, E., loc. cit., formulas (1.1), (3.3). 
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Allowing r to approach zero, from (6), (1/7) f2, ¢(x) dx = —1, so 
a ee | ~nba(*)h 
| + — P-) Tap = ~2e0V meat (1 — 7) 2 af) p. 
Therefore, from (6), 
= oa dx ae P 1 + N (3 be za)" 
| 5 (p. p-) 1—rz/b 2reV ba {t II. 1 — rl, i (9) 
Now let 


«Fs, 2 ES he Tt, 2 Be 6 


Then 


I]. = rt)/Q — rh) = 1 = (, — Ly) + 5°: — Ly - (1 - 14) + 00°), 


1 
so 
N 


Tl. G@ —7t)'70 — rly”? =1- 5 nT, = 


1 
+ 8 r((T, — L,)’ — °T, — L,)] + 06°). 
This may be inserted into (9), expanding the right side in powers of r. Expanding the 
left side in powers of r and equating coefficients, we have 
L = —xpV’*d(T, — La, 


(11) 


M= mpV'b'[(7; we L,)’ ~~ 2(T = L,) Ja, 


| 


for the aerodynamic lift and moment respectively. 

Suppose the slots are symmetrically arranged about the centerline, and suppose ¢ 
is the fraction of the total roadbed width occupied by lanes. Then T, = L,, T, — L, = 
2e, SO 


L = —2zpV’bea, M = xpV’b’ ea. (12) 


This proves the following theorem. 

Theorem. A flat plate airfoil slotted symmetrically about its centerline has the same 
steady-state aerodynamic lift and moment as an unslotted flat plate of chord equal to 
the total chord of the slotted plate minus the total width of the slots. 
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NOTE ON THE THERMAL STRESSES IN A LONG CIRCULAR 
CYLINDER OF m + 1 CONCENTRIC MATERIALS* 


By B. E. GATEWOOD (AF Institute of Technology) 


It has been shown that the stresses and displacements in a long cylindrical body‘ 
can be expressed in terms of two functions U and V which satisfy the equations 


V‘U = 0, V’V = kT, (1) 


where 7’ is the temperature change, and k = Eqga/(q — 1) with E Young’s modulus, 
q the reciprocal of Poisson’s ratio, and @ the linear coefficient of thermal expansion. 
The stresses are as follows 
o, = U,, — Vy; o, = U,, — Vaz; Te, = Vi, — U., 
(2) 
f., = 0, fT, = 0, o, = Ee, — kT + V’°U/q, 


where the constant strain e, is determined from the requirement that the total normal 
force on the cross-section of the cylinder be zero. The function U can be expressed’ in 
terms of two analytic functions ¢(z) and H(z) as 


U=2¢@+2¢2+H@+H@), (3) 


and the function V may be taken as a particular integral of V?V = kT since the solution 
of V’°V = 0 can be included in H(z). 

The boundary conditions for the case of a circular cylinder composed of m + 1 
different materials concentric to each other, each of which has different HE, g, and a, 
can be obtained by taking the displacements and normal stresses continuous on the 
junction surfaces and the normal stress zero on the outside surface (let S,, be the inside 
material, S,,, the adjoining material concentric to S, , and S, the outside material; 
let C, of radius rp be the outside boundary and C; of radius r; the boundary between 
S,; and S;_, with 7 = 1, 2, --+ , m): 


do(t) = t pi(t) + Yo(t) -_ Y Jo On Ge, (4) 
(1 + My)o;(t) + Nibis() + #10) + Hi 
(5) 
= R,t+[Y]; , on C; withj = 1, 2, --- ,m, 
M ¢;() + (1 + N 0; i(t) + t o/_,(t) + v;-1() 
(6) 
= Rt + [Y];-.,o0nC;, 
where z = tonC;, Y = V,+ iV,,G = Eq/2(q + 1), ¥(@) = A’ (2), and 
M; = 4G';_1(q; mas 1) qi(G; ‘= G;-1); J = l, 2, om 
—N; = 4G (q;-1 a 1) gi-i(G; = G;-1), (7) 


R; = 2e,.G,G; (Qi -1 — 9:)/ 9:9: -1(G; as G;-1). 


*Received Sept. 22, 1947. 
1B. E. Gatewood, Thermal stresses in long cylindrical bodies, Phil. Mag. (7) 32, 282-301 (1941). 
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Since the functions ¢; and y; are analytic, they can be represented by a Laurent 


series in S; (7 = 0,1, --- , m — 1) and a power series in S,, : 
:(2) = Do Gin + idjnle” + DS (Cin + tdjne”, 
n=O n=l 


(8) 


) 


yw (Cin + tf inde + > (Jin + th;,)z ", 


n=0 


¥;(z) 


With Cnn = Ann = Inn = hm, = 0. Also, Y; may be expressed in a ‘series of powers of 
z such that on C; (7 = 0, 1, --- , m) 


[Y]; = } (Ain + 1B;,)t" + Z (Cin + tD;,)t, 
; : (9) 
(Yli- = Do in + Fadl + DL Gin + tH dt. 


If these expressions be substituted in the boundary conditions (4), (5), and (6) and 
the resulting expressions together with their conjugates be integrated over the respective 
circles by use of the Cauchy integral formula, then a system of equations for determining 
the desired analytic functions is obtained. In some cases the functions can be deter- 
mined directly while in others they can be obtained by equating coefficients of like 
powers of z. For the general case the coefficients in Eq. (8) can be expressed in terms 
of the coefficients in Eq. (9) by means of recursion formulas. The coefficients ajo + tbjo , 
€io + tf, , bs, , and h;, may be neglected since they contribute nothing to the stresses. 
a;, and g;,; are determined by 


2a, + 9oi/Te = Ay, Im = 0, 
(2 + M;)a;; + N a; aat 9i1/7; - R; + Ain , J =1,°---,m, (10) 
Maj, + (2 + N;)aj;-1.1 + 9i-1,1/7; = R; + Ey, . 


Define Win = Gin + 1D in » Vin = Cin + id in > Yin = Cin + Uf in » Zin = Gin + thin , W in = 
A;, + 7B;, , ete., whence w;, and z;, are given by 


Wor + Zo2/To = Wor, ua = O, 
(1 + M,)wj. + Njwj-1.2 + Z2/7§ = Wis, j=l,-->,m, (11) 
M wy. + (1 + N,)j-1,.0 + 2)-1,2/75 = Via. 
The remaining coefficients w,, , £,,.-2 » Yi.n-2 , and z;, (n > 3) are given by (p; = r;4:/1r;) 
—M ja:p;"Wisrn +L + Ms — (L + ONissp" Join + Niwj-1.0 
— (0 — rE ;.n-0(1 — pi) = Win — PF Yisin» 
—M 5 :Fj41.»-2 + [CL + Mipi* — (Ll 4 Nios)Ein-a + Nive “Fir n-2 


2n-2/..-2 2n-4 F . 
+ nies (Di — Win = PF Xjn-2 — Zisre-a, Jw ++, ™, 
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Inn =O, Zan = O, Zon = 1o"|Won — Won + (n — 2) “"Xo.n-2], i 

Zon =n'[Y,, — Mw, — (1 + Ni)won + (n —- es SR | 

Zz, = r"(Win — (1 + My)wjn — Nywjeaan + (mn — Wri"; n-2)), 

Yoa-2 = ZX. 0 = 0 "Wes ~ Boca), 

Yo..-2 = r-?7Z, no — Mids.n-2 — (1 + Ni)Z0.n-2 — 02” "Worl; ' 

Yin-2 = an) — (1 + M;)%;.n-2 — Nidj-1n-2 — N7j"* Win). | 
For the case of the temperature a function of the radius the functions become very 


simple with the series in Eqs. (8) and (9) reducing to one term. With the temperature 
as 7T,(r) in S, , the particular integral of V’V = kT’ yields 


Il 


7. 


z. f rT’,(r) dr, 


2¥, =k [ rat) dr + hn [Tartar +++ +k; [| Par, 
“0 “rm “Cita 
whence, using Eq. (9), [Y]; = Lt = Ant = [Y];-1 = Zyt with r?L,; = [zY;],-,, . Further, 
¢,;(2) = a,,z and ¥;(z) = gj:/z, where a;, and g;, are given by Eq. (10). The function 
U is U; = aj’ + g;; log r and the radial and tangential stresses in polar coordinates 
are 
24;1 + (Gin = zY;) r’, T +0; = 0, 


O;; 


oo; = 2a;, — (gir _ 2zY,)/r’ _ k,T; ’ 


Eye, — k;T; + 40;:/q; . 


0; 


ON THIRD-ORDER CORRELATION AND VORTICITY IN 
ISOTROPIC TURBULENCE* 


By F. N. FRENKIEL (Cornell University) 
The fundamental equation of the propagation of the correlation function obtained ‘ 
Yr 2 , hd bd . . . ' 
by Th. von Kérm4n and L. Howarth’ for isotropic turbulence in an incompressible 
viscous fluid, is written 


a(u”f) ~ 12) an(ah 4 ) = er (=f 4 2s) 
ry + 2(u") ar os : h} = Qu ar “ a a (1) 


*Received Nov. 19, 1947. 
1Th. von Karman and L. Howarth. On the statistical theory of isotropic turbulence, Proc. Roy. Soc. 


London (A) 164, 192 (1938). 


ot 
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For small distances r, the longitudinal correlation function f may be represented 
by the relationt 
_ 1 2 ” -_ 4f1V 
and using the two relations 
\* = —f’(0) and A,” = av) 


where \ is Taylor’s dimension of the smallest eddies, there will be obtained for small r 


the following expression: 


fei-s+eo5 (2) 


Im _ 192 (3) 


and multiplying by \*/r’, one has, for small r, 


ad, 5r(,_® ‘er (2 4) - 
* 4 55(1-¥) +2 ne * tay dls (4) 


r 
This equation gives the law of the propagation of the dimension \. It may be also 


written in the following form 
dh, 4, _ e ) 1a vn 
rt," E ~~ atl wy (4a) 
and after integrating this differential equation we find 
5 (x? ) 1 a l (2) 
"™ E (s ~ 1) O87 at JNK 6) 


(u”)' /2 


Vv 





where N, = 


is the number of Reynolds of microturbulence. 

Let us consider the decay of isotropic turbulence and assume with von Karman 
the case of “large’’ Reynolds number of turbulence, i.e. the case when the correlation 
functions f(r, t) and h(r, ¢) are independent of the viscosity », except for small r. In 
addition, von Karman assumes that the correlation curves preserve their shape and 
that only a seale L changes. As a result the relation’ 

uy" 6) 

Ly ; 
is found, where A is a numerical constant and L a dimension representing the scale of 
turbulence 


+To avoid confusion between values at an initial moment é which appears later and values at the 
origin, we shall use in this paper f’’(0) for f” at r = 0 instead of fy’’, h’’”’(O) instead of ho’”’, and so on. 
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L. G. Loitsianskii’ shows, that when the rate of decrease of f, h and df/dr exceeds 
the rate of increase of r*, there exists the following theorem of the conservation of the 
disturbance moment: 

u® | s'f(@) ds = uw? L** = const. (7) 
“0 

As the disturbance length L* is given as a function of the fourth moment of f and 
as the shape of this curve is assumed to be preserved for large r, we may use L* in place 
of L, assuming that L*/L is constant. Thus the equation (6) becomes 


(u)/*y? Z 
L*y 


A, (6a) 


where A, is a new numerical constant. 
After eliminating L* between (6a) and (7) it is found that 


(uJ? = 


(8) 


o|— 


v 


with C a numerical constant. Eliminating now \°/v between the last relation and the 
equation (3) and integrating, we find 


(uy? — (us?) = 7C(t — &) (9) 


where u” = u/? at a chosen initial moment ¢ = ¢t, . C is found as a function of the initial 
values uj? and ), using equation (8), and as a result the following equation of decay is 
obtained: 


—— = yp —10/7 
u"?=u's [7 — (t — t) + | (10) 
It will also be found that 
iin Lal 7 ees i} (11) 
and 
= a 7 3 (t — to) + | (12) 
Xo 
or ? = r+ T(t — to) (12a) 


In another paper® these results are compared with experimental measurements of 
the decay of turbulence, and the agreement with equations (10) and (12) seems quite 


satisfactory. 
The variation of the Reynolds Number of microturbulence is easily found using 


equations (10) and (12) and is written 


—3/14 
Ny, = Mol 7 xz (¢— &) + | (13) 


2L. G. Loitsianskii. Some basic laws of isotropic turbulent flow, N. A. C. A., T. M. 1079, 1945 (Trans- 
lated from Rep. No. 440 of Central Aero-Hydrodynamical Institute, Moscow, 1939). 
8. N. Frenkiel. The decay of isotropic turbulence, (to be published). 
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For the decay of vorticity it will be found 
pa Gall y -17/7 
o 
with w? = 15u”/d", the value of the mean square of the vorticity, w. being its value at 


the initial moment f, . 
Differentiating equation (12a) it will be found 





1 dy’ 

-— 7 (15) 
which combined with equation (5) gives for the third order correlation coefficient the 
expression 

_|5_ 17] 1 (z) 
“< F *d z| Ny \n (16) 


In this relation the third order correlation coefficient h at the moment ¢ is given as a 


function of A, A,, and N, at the same moment t. 
Equation (16) represents a parabola of the third order to which the triple correlation 


curve h(r) is tangent at its origin. 





The last equation gives P 
aie Ce iz 6 
It should be recalled that 
? 7 f'"(0) 





i _ 15 "0)}° 


depends on the shape of the correlation curve. Using the notations of Batchelor and 


Townsend* S = —)*k’"(0); G = a*f'"(0) = f'"(0)/[f’O)? and remarking that for 
isotropic turbulence’ k = —2h and k’’(0) = —2h’’’(0) it will be found from (17) that 
fake 

G=t5™S (18) 


giving a relation between G, S and N, at the same moment t. 

Consider the case of decay of isotropic turbulence and assume that G = G, is con- 
stant during the decay. This assumption agrees with the results of the experiments of 
Batchelor and Townsend.* The constant value of G, depends upon the experimental 
conditions and may be different in each experiment, depending on the shape of the 
correlation curve. Thus during the decay 





CF 2 ( Yy st) 
S = N, Go 14 (19) 
and at the initial moment f, 
3/14 
S = s| 7 “3 (t — t) + | (20) 
0 


‘G. K. Batchelor and A. A. Townsend. Decay of vorticity in isotropic turbulence, Proc. Roy. Soc. 
London (A), 190, 534 (1947). 
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with S, = z Go - ay (21) 
N,.0 14 
The decay of the vorticity is given by the relation’ 
des” ~ ~72\3/2 w ¢ 
a 70h’""(0)-(u”)*’* — 10v 2 (22) 


and it will easily be found, using the relations given before, that 


dus? 7 4 ‘ G Eraser 

—— oe ere (wo, and 8, re 2 De. Nye (23) 
dt 3(15) : Nj.0 

In the decay of vorticity, the rate of production of vorticity by diffusive stretching of 

the vortex tubes and the rate of dissipation of vorticity by viscosity are each proportional 

to (w)**”"". The ratio of the production of vorticity to its dissipation is constant and 





equal to 








SN) — 51 1 (Q4 
; ~*~ fe ated 
The equations (23) and (21) give 
dex” _ 17 (ay ~2\24/17 
dt - (15)'” Ny. (w ) (25) 


which could be found directly by differentiating equation (14). 

There exist some differences between the conclusions drawn here for isotropic turbu- 
lence at “large’’ Reynolds number of turbulence and the conclusions of Batchelor and 
Townsend. In particular, in the present paper it is found that the decay of vorticity 
is proportional to (w”)**”"” and not to (w’)*”. The factor G is considered to be constant 
as in the Batchelor-Townsend paper but N, and S are variable. For the case studied 
here, G is given as function of S by the equation (18), which is of the same form as that 
found by Batchelor and Townsend but instead of 30/7 the constant is equal to 51/14. 

The experimental results of Batchelor and Townsend‘ seem to agree only very 
roughly with these relations. However, if account is taken of probable inaccuracies in 
the experimental data, such as the imperfect isotropy behind grids and errors due to 
finite lengths of hot-wires (especially as this concerns measurements of \), it appears 
that the agreement may be satisfactory after all. 


BOOK REVIEWS 


Relaxation methods in theoretical physics. By R. V. Southwell. Oxford, at the Clarendon 

Press, 1946. vii + 248 pp. $7.00. 

The relaxation method of approximate numerical solution of systems of equations was first intro- 
duced in 1936 by R. V. Southwell in connection with the solution of engineering problems arising in the 
field of structural design. In questions of the loading and deflection of complex structural frames, one is 
struck by the inherent precision of the conventional analytical methods of attack and the inherent lack of 
precision in the given data defining a problem. A number of people devised numerical schemes of calcula- 
tion by which the precision of the answer could be made as good as desired and thus made comparable 
with the given data. Several of these methods involved successive relaxation of constraints. To Southwell 





— 


Yawn 
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however, belongs the credit of realizing the great generality of this method of approach. In steps, he and 
his team of computers have solved more and more complex problems. 

In the book ‘‘Relaxation Methods in Engineering Science,”’ Southwell has presented in some detail, 
problems involving systems of algebraic equations and ordinary differential equations, including eigen- 
value problems, which arise in engineering structures and vibrations problems. 

In the present volume, the treatment of problems by the relaxation method is extended to include 


partial differential equations from various branches of classical physics. All of the problems are included 


in the general form 


0 ( Og 0 
—({y—]} = Z. 1 
Ox 3) Oy \" dy () 


where 
¥ = V(2, y, >, b: » by), Z = Z(t, Y, o, bz » by)- (2) 


The first chapter is devoted to a discussion of general ideas and formulas. The relaxation method is 
explained and the finite difference expressions for differentials of various orders are derived. Partial dif- 
ferentials are then treated, and in particular the finite difference approximation to the Laplace operator 
is derived for triangular, square and hexagonal nets. Numerical integration and interpolation are also 
briefly treated in this chapter. 

The remainder of the book, chapters II through VI are devoted to successively more difficult cases of 
(1) and (2). The added difficulties are of two classes: (a) more difficult boundary conditions, and (b) more 
complex forms of (1) and (2). 

In chapter II, Equation (1) with y = 1, Z = Z(z, y) is discussed in detail. (No significant numerical 
simplification results if Z = 0.) The relaxation pattern, the basic idea of relaxation, block relaxation for 
speedier approach to the solution, advance to a finer net, and the physical significance of these techniques 
n terms of the “tensioned net’’ are all carefully presented. The chapter concludes with several torsion 
problems, some of which are left for the reader to try. 

In chapter II, the boundary conditions are specified values of the wanted function on specified 
rectangular boundaries. In chapter III, the treatment is broadened to include specified normal gradients 
and irregular boundaries. Several problems with known analytical solutions are solved to illustrate the 
accuracy attained with relatively course nets and then several much more difficult problems involving 
the torsion of irregular shafts and magnetic fields in electric machines are solved to illustrate the power of 
the method. 

Chapter IV deals with general questions of conformal transformation. The numerical determination 
of the real and imaginary parts of the complex function required to transform any given domain into a 
circle or rectangle is considered. A treatment for suitable removal of the logarithmic singularity at the 
center of the circle is given as is the determination of the complex conjugate of a function. It is shown 
that the most accurate determination of a complex conjugate generally involves a separate relaxation 
for each. 

Cases of (1) with y = ¥(z, y), of which the Laplace operator in polar coordinates (7, z) is a typical ex- 
ample, are treated by the use of relaxation patterns which vary from point to point. Chapter V presents 
the necessary theory and gives a considerable number of illustrative problems from the fields of elasticity, 
heat conduction, and lubrication. This chapter concludes with two problems in which y is a function of @ 
and its derivatives as well as x, y. A lubrication problem with variable viscosity and the flow of a gas 
through a nozzle are the illustrative problems treated in detail. 

In the final chapter (VI) “Problems Involving Boundaries or Interface Not Initially Known” are 
treated. Usually boundary conditions consist of two parts, a given boundary and a given relation between 
the function and its normal derivative along this boundary. Instead, one may choose two other suitable 
conditions say the value of the function and the value of its normal derivative but omitting the shape 
of the boundary itself. Such problems arise in a natural way in many physical problems. The problems 
discussed involving plastic torsion, percolation of water through soils, the waterfall, and the wake behind 
a cylinder are all cases in which analytical methods have given few or no solutions. 

No exaggerated claims are made for the ease by which the problems can be solved by the relaxation 
method. It is in fact pointed out that problems like those of chapters V and VI especially involve consider- 


able labor of a type which does not appeal to many persons. It is correctly noted however that problems 
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of these and even more difficult types (not yet released from wartime publication restrictions) can be 
solved by the relaxation method, whereas they are solvable in no other way at present. 

Throughout the book, it is emphasized that the only way to appreciate the significance of the relaxa- 
tion method of numerical solution of partial differential equations of mathematical physics is to do 
problems. Thus each chapter presents the physical theory and finally accepted answer to a number of 
problems, leaving the numerical work as exercise for the reader. 

Numerous tables of finite difference formulas for differentiation and integration are given in 10 pages 


of appendix. 
ot appe: Howarp W. Emmons 


Conduction of heat in solids. By H. 8. Carslaw and J. C. Jaeger. Oxford at the Clarendon 

Press, 1947. viii + 386 pp. $8.00. 

This book is a very extensively revised and enlarged version of the senior author’s “Introduction to 
the Mathematical Theory of the Conduction of Heat in Solids” which appeared in 1921. It treats very 
thoroughly the best available mathematical techniques for the analysis of heat conduction and includes 
a large number of solutions to specific problems. 

The first chapter presents a discussion of the basic partial differential equation, and the mathematical 
formulation of initial and boundary conditions for typical systems. In the following eight chapters the 
Fourier methods of analysis are applied to determine transient temperature distributions under various 
conditions in infinite and semi-infinite solids, slabs, rods, rectangles, rectangular parallelepipeds, circular 
cylinders, wedges, spheres and cones. In Chapter X, the method of sources and sinks and the method of 
images are developed and applied to problems not easily solvable by the methods of the preceding chap- 
ters. Much, though by no means all of the material in these ten chapters was contained in the earlier 
book; the remainder of the book is largely new. 

The next four chapters present a clearly written discussion of the Laplace transform and a thorough 
treatment of its application to the solution of heat conduction problems. This very powerful method, 
which was developed after the appearance of the earlier book, replaces the “method of contour integrals” 
discussed therein, and furnishes solutions to many complex problems which do not yield readily to the 
more elementary techniques. Many of the original results obtained by the authors are reproduced or 
presented for the first time in these chapters. The Laplace transform is also used to obtain the Green’s 
functions for heat conduction in solids of various typical geometries. A final chapter includes for com- 
pleteness a very brief treatment of steady state heat conduction. 

Many graphs and tables not contained in the earlier book and useful in making calculations for 
engineering applications are included for ready reference. A very brief discussion of the available numer- 
ical and graphical methods is also given. Many references to the literature are made, but there has been 
no attempt to compile an exhaustive bibliography. The theory developed is a linear one, i.e. no attempt 
has been made to treat the problem of variable thermal properties. A minor defect is that the table of 
contents and the index are possibly over-concise, so that a page-by-page search is sometimes necessary 
to locate a desired result. The book as a whole is a clearly written, complete and excellent treatise on the 


subject. , . 
q S. A. ScHAAF 


Tables of supersonic flow around yawing cones. By the staff of the Computing Section, 
Center of Analysis, Massachusetts Institute of Technology. Under the direction of 
Zdenek Kopal. Cambridge, Massachusetts, 1947. xviii + 321 pp. $5.00. 

In the introduction to this volume, the theory governing the supersonic flow of a gas past a slightly 
yawed cone is developed. In this development, the Taylor-Maccoll solution for a symmetrically oriented 
cone forms the zero order terms of a perturbation expansion with regard to the yaw angle e. The equations 
relating to the terms of order ¢ are derived. The original detailed theory which has not been previously 
published is said to be forthcoming in a periodical. 

The problem reduces to the determination of the coefficient of « in the perturbation development 
and the body of the book consists of a compilation of the numerically obtained results. The range of cone 
angles is 5° < 3, < 50°. As in volume I (the corresponding tables for the unyawed cone) results are 


obtained for a large number of wave angles corresponding to each cone angle. 
G. F. Carrier 
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Probléme général de la stabilité du mouvement. By M. A. Liapounoff. Princeton University 
Press, Princeton, 1947. 273 pp. $3.50. 


Number seventeen in Annals of Mathematical Studies, this book is a translation from the Russian. 
Published in Russian in 1892 by the Mathematical Society of Kharkow, the French translation appeared 
in Annales de la Faculté des Sciences de Toulouse, 2nd ser., vol. 9 (1907). 

The book is divided into three chapters, the first of which is labeled Preliminary Analysis; it includes 
definition and discussion of the general problem of stability. Chapter I also contains a discussion of sys- 
tems of linear differential equations and differential equations of the disturbing motion. Chapter II is 
entitled Study of Steady Motion and it includes more specific treatment of the differential equations of 
the disturbing motion for various cases of steady motion. Chapter III is concerned with the Study of 
Periodic Motion and includes analysis of the characteristic equation as well as the differential equations 
of the disturbing motion. 

This work encompasses the classical treatment of stability for continuous motion of disturbance. 


There is not any discussion of dynamical processes subject to discontinuous impulses. 
; R. TRUELL 


Theory and application of Mathieu functions. By N. W. McLachlan. London: At the 
Clarendon Press, 1947. xii + 401 pp. $12.50. 


In the first two-thirds of this treatise the developments of the solutions of the Mathieu equation are 
carried out and various techniques for the evaluation of the necessary parameters are given in some detail. 
These techniques, both numerical and analytic, are amply illustrated by appropriate numerical examples, 
and their ranges of validity are indicated. The solutions are developed in the form e“*¢(z) where ¢ is found 
as a trigonometric series, in the form of Maclaurin series, and in series of Bessel functions. The eigenvalue 
problem associated with the Mathieu equation is treated. Orthogonality conditions and several integral 
identities are established. The asymptotic representations and zeros of Mathieu functions are discussed. 
The Hill equation is briefly treated. 

In the second part, various physical problems are solved in terms of Mathieu functions and the 
results are interpreted. Loud speaker distortion, the oscillations of an elastic column under periodic end 
thrust, the vibration of a membrane with an elliptical boundary, the viscous flow past an elliptic cylinder, 
and other problems, are used to illustrate the large class of problems in which Mathieu functions arise. 


G. F. Carrier 


Tables of spherical Bessel functions. Prepared by the Mathematical Tables Project, 


National Bureau of Standards. Vol. II. Columbia University Press, New York, 1947. 

xx + 328 pp. $7.50. 

In this continuation of Vol. I of the same work, the function (#/2zx)/*J (x) is tabulated for the half- 
odd-integer values 29/2 < | »| < 61/2. In the range 29/2 < | »| < 43/2, the argument z ranges from 0 
to 10 in steps of .01 and from 10 to 25 in steps of .10. For 45/2 < | »| < 61/2 the range of z is 10 to 25 
in steps of .10. For } < » < 61/2 and for —29/2 > vy > —61/2, the function 2’T(» — 1)J(x)/z” is tabu- 
lated for 0 < x < 25 at increments of .1. The number of significant figures is at least seven. Finally, 


“ 


some of the zeros, maxima, and minima, of the half order Bessel functions are tabulated for —29/2 < 


vy < 45/2 
G. F. Carrier 


Vector and tensor analysis. By Louis Brand. John Wiley & Sons, Inc., New York, 1947. 
xvi + 439 pp. $5.50. 

This text contains a very good introduction to the classical theory of vectors, motors, dyads, tensors, 
and quaternions and some of their applications to geometry and physics. Throughout the text, the main 
emphasis is on developing the student’s manipulative ability for dealing with these powerful tools. This 
is accomplished by offering very complete explanations of techniques and using numerous illustrative 


problems throughout the text. 
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The first section deals with Gibbs’ Vector Analysis (Chapters I, III, V, VI). The dot, cross, triple 
sealer, etc. products and reciprocal basis are introduced in Chapter I. As applications, the author offers 
proofs of various theorems in classical plane projective geometry and vector interpretations of homogene- 
ous, line, and plane coordinates. Following the introduction of vector differentiation and integration 
(Chapter III), the author treats the theory of curves (the Frenet formulas and applications). As a further 
illustration of vectorial methods over one dimensional fields, the kinematics of a particle is discussed 
The portion of the text which treats vector fields of two and three dimensions (Chapters V, VI) discusses 
such topies as the gradient, the curl, the divergence, vector derivative identities, curvilinear coordinates, 
surface vectors, the Green and Stokes theorems and some elementary potential theory. It is the reviewer’s 
opinion that the author uses the cross product in several instances where use of the dot product would be 
more advantageous. Thus, in deriving the Frenet formulas, the author works with the cross product. 
Here, the dot product furnishes an easier derivation and one that can be generalized to n-dimensional 
manifolds. 

The algebra of motors is studied in Chapter II. This work is used in only one application (Chapter 
III, pp. 117). Thus, this topic may be included or excluded in any course based on this text without af- 
fecting the continuity of the remaining material. — 

Tensors and their applications form the basis of Chapters IV, VII, VIII, IX. The classical Gibbs’ 
theory of dyadic algebra is discussed in Chapter IV. Topics such as conjugate dyadics, decomposition of 
dyadics, symmetry, invariants of dyadics, relation between dyadics and matrices are treated. Some use of 
these ideas is made in Chapters V and VI. In the main, the use of dyadics in these last mentioned chapters 
serves to generalize the Green and Stokes vector theorems. However, fairly extensive use of dyadics is 
made in the treatment of Hydrodynamics (Chapter VII) and Surface Geometry (Chapter VIII). The 
section on Hydrodynamics discusses the Euler and Lagrange equations, the Helmholtz equations, and 
the Kutta-Joukowsky formula. In the Chapter on Surface Geometry, the author treats curves on a sur- 
face by means of the Darboux moving trihedral. Geodesic curvature, geodesic torsion, normal curvature, 
lines of curvature, asymptotic lines are among the many topics introduced. Further, the first and second 
fundamental forms, geodesics, the mean and total curvature, the Gauss-Bonnet formula, and the Gauss- 
Codazzi relations are all discussed in some detail. Chapter [IX is concerned with Tensor Analysis. Trans- 
formation laws of vectors and tensors, multiplication and contraction of tensors, tensor densities, the 
permutation symbol, and the generalized Kronecker tensors are examined. By expanding geometric 
objects (dyads, triads, etc.) in terms of their tensor components and a system of base vectors and by 
determining the differentiation law of the base vectors, the author introduces covariant differentiation 
The Riemann-Christoffel and Ricci tensors are discussed. The reviewer would like to make two criticisms 
of the work in this section. First, the.use of dyads in the section on Surface Geometry (Chapter VIII) 
needlessly increases the amount of technique which students must develop and, further, detracts from the 
unity of the text. The same results can be obtained by the proper use of tensors (components of the 
dyads). Secondly, in deriving the Kutta-Joukowsky formula, the author passes from vectors to complex 
numbers by a vague equivalence scheme. It is known (see A. Wundheiler: Are Complex Numbers Vectors?: 
Bull. of Am. Math. Soc., Vol. 46, No. 1, Jan. 1940, pp. 57) that “complex numbers are not isomorphic 
with vectors but are vectors in a metric plane with a distinguished direction.’’ Thus, it appears desirable 
to show clearly the relation between complex numbers and vectors or else to treat complex numbers 
independently of vectors. 

In the final Chapter X of the text, Quaternion Algebra is developed. Some applications are given to 
rigid body motions. 

The author has presented in an interesting and readable manner several methods of treating geo- 
metric quantities. As a result, this book should form a valuable addition to the existing texts on Vector 


and Tensor methods 
N. Cospurn 


Sequential analysis. By Abraham Wald. John Wiley & Sons, Inc., New York, 1947. 

xii + 212 pp. $4.00. 

One of the most notable recent developments in statistical theory is sequential analysis, a sampling 
procedure in which the size of the sample is not fixed in advance, but depends on the results of the succes- 
sive observations. The basic results are due to Professor Wald, and have been published in two papers 
in the Annals of Mathematical Statistics; this book is essentially an expanded and somewhat reorganized 
version of Professor Wald’s papers, together with some additional material. 
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The book consists of three parts and an appendix. Part I contains background material and most of 
the theory, inciuding a discussion of the sequential probability ratio test, the most important sequential 
procedure. Part II is a detailed description of sequential procedures for testing certain hypotheses of 
practical importance; the hypotheses treated concern the parameters of binomial and normal populations. 
Part III is a brief discussion of some possible sequential approaches to the problem of multi-valued de- 
cisions and estimation. The appendix contains those parts of the theory which require some advanced 
mathematics; the rest of the book is not beyond the understanding of one acquainted with elementary 
calculus. 

Two features of the sequential method deserve especial mention. 1) The sequential procedure re- 
quires on the average fewer observations than comparable fixed-sample-size procedures, the saving in 
some cases of practical importance being as much as 50%. 2) Sequential theory is much simpler than that 
for samples of fixed size. The book is a clear presentation of an important topic, and should be invaluable 


to practical statisticians. 
Davin BLACKWELL 


Guide to the literature of mathematics and physics including related work on engineering 
science. By Nathan Grier Parke III. McGraw-Hill Book Co., Inc., New York and 
London, 1947. xv + 205 pp. $5.00. 


This ‘‘Guide”’ is planned for the “large group of people who cannot maintain an independent aware- 
ness of the mathematical and physical reference literature but who can make effective use of a classified 
guide. This group includes scientists, engineers, librarians, and students . . . Often the techniques that 
overcome obstacles in one field can be effectively used in another. For example, linear electrical networks 
are the most intensively studied vibration systems. By means of electromechanical analogues, all these 
results from the electrical field are applicable to the mechanical field.’’ The bibliography of approximately 
1800 entries, classified under some 150 headings, is preceded by over 70 pages of detailed discussion of 
the principles of reading and study, self-directed education, literature search, and too little on the use of 
periodicals. Suggestions are offered on problems of library technique and use. 

As the intended audience is so inhomogeneous, both in interests and in level of accomplishment, the 
author has necessarily assumed almost complete lack of knowledge on the part of the reader. An experi- 
enced worker, however, would do well to read all of the, what might appear to him as a somewhat 
pedantic and overdetailed, discussion for the valuable hints it contains. Unfortunately a difficulty arises 
in the use of the bibliography whose primary functions should be to tell the reader about the essential 
treatises in each field, books which bring the information up to date, and also periodicals in which reviews 
and research reports are published. In some cases, the author has done this carefully, but in many classi- 
fications, the research worker or engineer who is branching out, possibly temporarily, into an unfamiliar 
field is confronted with too many books about which too little is said for him to make an intelligent 
selection 

In its present form, the “Guide” is most useful to those who have an elementary understanding of 
the topics considered and who require detailed or advanced knowledge. It is to be hoped that frequent 
revisions or supplements will be issued to keep the “Guide” valuable. As an example, illustrating this 
necessity, Applied Mechanics Reviews will soon be published as a periodical covering the current world 
literature in many of the classifications. Obviously, no mention of this valuable addition to the review 


journals can appear in the present “Guide.” 
D. C. DruckER 


Table of the Bessel functions J)(z) and J ,(z) for complex arguments. Prepared by the Mathe- 
matical Tables Project, National Bureau of Standards. Second edition. Columbia 
University Press, New York, 1947. xliv + 403 pp. $7.50. 

The preface states that since the publication of the first edition of these tables [reviewed in this 

Quarterly, 2, 276 (1944)] no errors have been reported in the tabular material. Accordingly, this material 

ould be reproduced from the negatives which were used for the first edition. Some corrections and revi- 

sions were made in the Introduction, however, which now contains the relations between the tabulated 


functions on the 45° ray and the ber and bei functions. 
W. PRAGER 
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Nomograms of complex hyperbolic functions. By Jérgen Rybner. Jul. Gjellerups Forlag, 
Copenhagen, 1947. xxxvi + 54 pp. Price Danish Kroner 24.00. 


This table of nomograms includes an introduction containing a note on the history of real and com- 
plex hyperbolic functions and notation, earlier tables and graphical representations of complex hyperbolic 
functions, the use of nomograms for: cosh (b + ja) = p + jg, sinh (b + ja) = p + jq, tanh (b + ja) = 
r | =rel?zt+tijy=r | 6 =rel*andR | w =1+r | as well as the use of nomograms for reflection 
losses and reflection phase shift and for resonance circuits and filters. Approximately ten pages of formu- 
las for circular and hyperbolic functions.include real functions, gudermannijan angle, integrals, differen- 
tials, series developments, functions of imaginary argument, functions of complex argument and inverse 
complex hyperbolic functions. 

A table of multiples of x/2 from 1 to 100 is present as well as some formulas for four terminal net- 
works and transmission lines. 

All of the material mentioned above is contained in thirty five pages (84 x 12). In addition there are 
fifty pages of nomograms as follows: 

cosh (b + ja) = p + jg 13 nomograms covering the range b = 0 to b = 4.0 nepers 
sinh (b + ja) = p + jq 13 nomograms covering the range b = 0 to b = 4.0 nepers 
tanh (b + ja) =r | @ 16 nomograms covering the range b = 0 to 2.0 nepers 


Eight additional nomograms of elementary functions related chiefly with transmission line and network 
calculations. 

This book of nomograms has been prepared with unusual care and accuracy, and its contents should 
be of considerable value to electrical and acoustical engineers and possibly to physicists. Any extensive 
numerical calculations involving impedances or reflection and transmission coefficients can be facilitated 


by means of these nomograms. 
Roun TRUELL 


























